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*5 ' Abstract 

This paper reports the Feynman-parametric representation of the vacuum-polarization function 

D : 

i-C ■ consisting of 105 Feynman diagrams of the eighth order, and its contribution to the gauge-invariant 

| set called Set I(i) of the tenth-order lepton anomalous magnetic moment. Numerical evaluation of 

>'.... 

' this set is carried out using FORTRAN codes generated by an automatic code generation system 

■ GENCODEVP./V developed specifically for this purpose. The contribution of diagrams containing 

in ; 

, electron loop to the electron g—2 is 0.017 47 (11) (a/7r) 5 . The contribution of diagrams containing 

| muon loop is 0.000 001 67 (3) (a/V). The contribution of tau-lepton loop is negligible at present. 

The sum of all these terms is 0.017 47 (11) (a/V) 5 . The contribution of diagrams containing 
^ . electron loop to the muon g—2 is 0.087 1 (59) (a/vr) . This is to be compared with the unpublished 

asymptotic analytic result (0.252 37 + 0(m e /m fl ))(a/ir) 5 . The contribution of tau-lepton loop 

to is 0.000 237 (1) (a/-7r) 5 . The total contribution to a^, the sum of these terms and the 

mass- independent term, is 0.104 8 (59) (a /vr) 5 . 

PACS numbers: 13.40.Em, 14.60.Ef, 12.39.Fe, 12.40.Vv 
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I. INTRODUCTION 



The anomalous magnetic moment g—2 of the electron has played the central role in testing 
the validity of quantum electrodynamics (QED) as well as the standard model. The latest 
measurement of a e = (g—2)/2 by the Harvard group has reached the precision of 0.24 x 10 -9 



[i 



l2|: 



a e (HV08) = 1 159 652 180.73 (0.28) x 10~ 12 [0.24ppb] . (1) 
: theory consists of QED corrections of up to the eighth order 



At present the best prediction o 

nil nn 

[3H5J, and hadronic corrections [6Hl2| and electro- weak corrections [13Hl5l| scaled down from 
their contributions to the muon g — 2. To compare the theoretical prediction with the 
experiment (JTJ) , we also need the value of the fine structure constant a determined by a 
method independent of g—2 . The best value of such an a has been obtained recently from 
the measurement of h/m-^, the ratio of the Planck constant and the mass of Rb atom, 



combined with the very precisely known Rydberg constant and m^K/m, |16l|: 

or -1 (RblO) = 137.035 999 037 (91) [0.66ppb]. (2) 

With this a the theoretical prediction of a e becomes 

a e (theory) = 1 159 652 181.13 (0.11)(0.37)(0.77) x 10~ 12 , (3) 

where the first, second, and third uncertainties come from the calculated eighth-order QED 
term, the tenth-order estimate, and the fine structure constant (j2J), respectively. The theory 
([3]) is thus in good agreement with the experiment ([1]): 

a e (HV08) - a e (theory) = -0.40 (0.88) x 10" 12 , (4) 

proving that QED (standard model) is in good shape even at this very high precision. 

An alternative test of QED is to compare the a of fl2]) with the value of a determined 
from the experiment and theory of g — 2 : 

a -\a e 08) = 137.035 999 085 (12)(37)(33) [0.37ppb], (5) 

where the first, second, and third uncertainties come from the eighth-order QED term, 
the tenth-order estimate, and the measurement of a e (HV08), respectively. Although the 



uncertainty of a _1 (a e 08) in (J5J) is a factor 2 smaller than a: _1 (RblO), it is not a firm factor 
since it depends on the estimate of the tenth-order term, which is only a crude guess 
For a more stringent test of QED, it is obviously necessary to calculate the actual value of 
the tenth-order term. In anticipating of this challenge we launched a systematic program 
several years ago to evaluate the complete tenth-order term 
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20|. 



The tenth-order QED contribution to the anomalous magnetic moment of the electron 
can be written as 



,(io) 



A^ 0) + ^ 10) (m e /m M ) + A { 2°\m e /m T ) + A^ 0) (m e /m^m e /m 7 



(6) 



where the electron-muon mass ratio m e /rria = 4.836 331 71 (12) x 10~ 3 and the electron-tau 
mass ratio m e /m T = 2.875 64 (47) xlO -4 17|. The contribution to the mass-independent The 
contribution to the mass-independent term A^ may be classified into six gauge-invariant 
sets, further divided into 32 gauge-invariant subsets depending on the nature of closed 
lepton loop subdiagrams. Thus far, 24 gauge-invariant subsets which consist of 2785 vertex 
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2lM24| . Throughout the paper the overall 



diagrams, have been evaluated and published 
factor (a /it) 5 is omitted for simplicity. 

In this paper we report the value of A^ contributed by a subset, called Set I(i), which 
consists of 105 Feynman diagrams obtained by insertion of proper eighth-order vacuum- 
polarization diagrams in the second-order anomalous magnetic moment M?,. These diagrams 
can be represented by 39 independent integrals taking account of various symmetry proper- 
ties. The evaluation of these integrals would be straightforward if the spectral function of the 
eighth-order vacuum-polarization were known. Unfortunately, it is not available at present. 
Thus we follow an alternative approach of expressing the eighth-order vacuum-polarization 
function U.^ s '(q 2 ) as a set of Feynman-parametric integrals and inserting them in the virtual 
photon line of the second-order anomalous magnetic moment M2 25]. 

Construction of the Feynman-parametric integral of the vacuum-polarization function 
and removal of subdiagram ultraviolet (UV) divergences by i^-operation 25| are described in 
Sec. [Til This scheme is implemented by an automated code generation system gencodevpTV 
developed specifically for this purpose. Incorporation of II^ 8 ^ in M2 is carried out in Sec. IHIl 
Since the .fT-operation subtracts only the UV-divergent part of the renormalization constant, 
additional removal of UV-finite parts of renormalization constants must be carried out to 
obtain the standard on-the-mass-shell renormalization. This is shown explicitly in Sec. HVl 
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(a)Type f 



(b)Type g 



(c)Type h 





v/W 



(d)Type i (e)Type j 

FIG. 1: Five types of diagrams that contribute to the eighth-order vacuum-polarization. 

Numerical evaluation of M^p is described in Sec. IVj where the first e in the superscript 
(ee) refers to the open electron line and the second e refers to the closed electron loop. The 
contributions of the muon loop and tau-lepton loop to the electron g—2, namely and 
Mgpg, are described in Sec. IVII The contribution of the Set I(i) diagrams to the muon g—2 



is described in Sec. IVHl Section IVIIII is devoted to the summary and discussion of this 
work. Especially, our result of the electron- loop contribution to the muon g—2 is compared 
to the prediction based on the renormalization group [26| and to the result obtained by the 



analytic-asymtotic expansion 27|, [28 ] . 



Appendix [A] describes the construction of Feynman-parametric integrals for M 2t p* . Ap- 
pendix [B] describes the on-shell renormalization scheme for the vacuum-polarization function. 
Appendix O gives intermediate renormalization of individual diagrams by the .fT-operation. 
Appendix [D] gives the divergence structure of quantities of sixth or lower orders. 

II. PARAMETRIC INTEGRAL OF VACUUM-POLARIZATION FUNCTION 



Diagrams that contribute to the eighth-order vacuum-polarization can be classified into 
five types according to their structures (See Fig. [I]). Contributions from the diagrams of 



Types 
Refs. I 



g and h, and j to the tenth-order lepton g — 2 have been evaluated previously in 



21. 



221 ]. respectively. In this paper we focus our attention on the remaining Type i, a 



set of 105 proper eighth-order vacuum-polarization diagrams, which is the most complicated 
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FIG. 2: ab(acbd) (cede): An example of sequential representation of a diagram of Set 1(g). 
one of diagrams shown in Fig. [I] and evaluate its tenth-order contribution Set I(i) to g — 2. 



A. Diagram representation 

In order to deal with diagrams which contain closed lepton loops as well as open lepton 
paths, we have to generalize the rules for the diagrams without closed lepton loop described 



in Ref. [19]. 



We begin by representing a diagram in terms of a sequence of symbols that characterize 
the photon lines by the following rules: 

1) Assign indices to photon lines, e.g. by lower-case alphabets, 'a', 'b', .... 

2) Identify a vertex by the index of photon line that is attached to the vertex. 

3) Read the indices of vertices along a lepton path (or loop) in a certain direction. (We 
adopt the reverse of the direction of the lepton propagator.) 

4) Enclose the sequences of indices of closed lepton loops by parentheses (but not indices of 
open lepton lines). 

For example, the tenth-order diagram with two lepton loops shown in Fig. |2] (which 
belongs to Set 1(g)) may be represented by a sequence, "ab(acbd) (cede)". 

This representation is not unique because there are several possible choices of assignment 
of photon line indices, cyclic permutations of vertices along lepton loops, and permutations of 
lepton loops and paths. To reduce the ambiguity we adopt the convention: the sequence for 
the open lepton path comes first, followed by the lexicographical sequences of closed lepton 
loops. The sequence within a loop is chosen also lexicographically, e.g. the sequence (dacb) 
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is rotated into (acbd). The photon line indices are taken from 'a' in order of appearance in 
the sequence. 1 

For diagrams describing a vacuum-polarization loop, which is our main concern, we adopt 
an additional rule that the two photon lines external to the vacuum-polarization loop are 
labeled by 's' and 't', whose Lorentz indices are /x and v, respectively. We also assume that 
the external momentum q flows in from the photon line £ t' (u) and leaves from the photon 
line 's' (//). The sequence of lepton lines in the loop is chosen to start from the index 's'. 

B. Algorithm to generate a proper lepton loop diagram 

We now present an algorithm for generating proper lepton loops of 2n-th order. A diagram 
of this type has a single lepton loop that consists of 2n vertices, 2n lepton lines, two external 
photon lines, and (n — 1) internal photon lines attached to the lepton loop. All lepton lines 
are directed, and two external photon lines are distinguished. 

The algorithm is as follows: 

1) A vertex to which an external photon line labeled by 's' is attached is chosen as the first 
element of the sequence. Assign the index '0' to this vertex, and assign numeric indices 
to other vertices sequentially along the loop in a certain direction. 

2) Another vertex is chosen to which the other external photon line labeled by 't' is attached. 
There are (2n — 1) choices of vertices. 

3) The remaining (2n — 2) vertices are made into (n — 1) pairs. Each pair corresponds to 
an internal photon line that connects the two vertices of that pair. There are (2n — 3)!! 
ways to construct (n — 1) pairs. 

Therefore, the total number of diagrams is (2n — 1) x (2n — 3)!!. 

Taking into account the time-reversal symmetry of QED and the symmetry by the per- 
mutation of Lorentz indices of external photon lines, we identify the equivalent diagrams 
with respect to the reversal of sequences and exchange of symbols 's' and 't'. In this fash- 
ion we obtain a complete set of topologically distinct diagrams with an appropriate weight 

1 This convention may still have ambiguities. However, it works for the diagrams with a single lepton loop 
that are discussed in the present article. 
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factor of the symmetry. In the case of the Set I(i) we have 39 distinct diagrams, which are 
shown in Fig. [3j 



C. Photon self-energy amplitude 



The momentum representation of the 2nth-order vacuum-polarization diagram G has the 
form given by the Feynman-Dyson rule: 



in^(g) = (-l)(-ze) 2 " J 



(2vr) 4 
Tr 



(27T) 



r 



m 



-r 



•7 



H+l 



m 



"7a 



n 



(7) 



where m is the rest mass of loop leptons, pi is the momentum flowing on the lepton line i, 
and kj is the momentum flowing on the photon line j. These momenta are given as linear 
combinations of the loop momenta Zi, . . . , l n , and the external momentum q. 

As a convention, the flow of the external momentum q is chosen as shown in Fig. HJ where 
each fraction q/2 flows on the upper- (lower-) semicircle that consists of lepton lines 1, . . . , i 
(i + 1, . . . , 2n), respectively. The function ([7]) is quadratically divergent, and we assume that 
the above expression is appropriately regularized by the Pauli-Villars regularization. 

We adopt here an approach that exploits gauge invariance of the sum IP" = ^ G Uq' 
which allows us to ignore the gauge- dependent part of ([7]). The gauge invariance ensures 
the identity 



q a U m {q) = 0, 



(8) 



which holds for Pauli-Villars regularized function IP". Differentiating it with respect to q^ 
we obtain 

dU au 



IF"(g) = -q 



(9) 



Since q^ is the external momentum, the order of g M and the integration over the loop momenta 
can be interchanged, as far as the integral is properly regularized. Thus we can write 

d^l\ d^l n 



m&b) = -(-i)He)% 



(2tt) 4 (2tt)^ 



_d_ 

dq^ 



■Tr 



7 



*>\ — m 



•••7 



h+i - m 



n 



• (io) 



PV 
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FIG. 4: Flow of extra momentum q. 



Carrying out the differentiation with respect to q^, we obtain 



'ng"fe) = -(-lX-ie) 2 ",: 2 " 



>i — m 



m 



(2vr) 4 (2tt) / 



2 J — m 



• • -7 



i+i - ™ 



n 



j v 3 



PV 



:n) 



where the minus (plus) sign in (=P7 M /2) is taken when the line j belongs to the upper (lower) 
semicircle of the diagram. 



The calculation can be simplified using the identity 25 ] 



+ mi)Y{i>i + m 

(p? - mfj 1 



2D? (pi + rm) 



(12) 



where D? is defined by 



i roc q 

* Jm 2 oq 



(13) 



We can now move the trace operation outside the k integration. With the help of Feynman 
parameters Zi associated with the line i, the denominators can be combined into one, and 
then the momentum integration can be carried out analytically. Now we bring back D- 
operators inside the z integral and obtain 

n-(g) = - (-£)" / ^ (£ ^) Tr \m +m)...T»...]±, (14) 
J j 

where 



(dz) = 5(1 — Zi) Y\ dzi, i = 1, 2, . . . , 3n — 1. 



(15) 
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U is a homogeneous function of z's determined from the topology of the diagram and V is 
defined by 



V = V - q 2 G, 



V = ^ Z i m2 

all lepton lines 



G= Z ^ 



(16) 



where the path V(fi, v) is arbitrarily taken to run between two external photon lines. In the 
present case, we choose the convention that the path V runs on the upper half part of the 
loop of diagrams shown in Fig. |H 

From Lorentz invariance and gauge invariance we have the general structure 

n GSk<?) = (?a»9" ~ I 2 9^) n§ n) (g 2 ) + (gauge-dependent terms). (17) 

The Lorentz scalar IF 2 ™) has only a logarithmic divergence so that the number of auxiliary 
masses in Pauli-Villars regularization can be reduced to one. The renormalization of subdi- 
agram UV divergences are carried out independently of each other. Gauge-dependent terms 
cancel out when all diagrams of Set I(i) are combined. The charge renormalization can be 
achieved by 



ng n V) = nr(9 2 )-nr(o). 



(2n)/ 2> 



- r (2n), 



(18) 



For the case n = 4, which is our concern, 11^ (g 2 ) can be expressed in the form 



(dz) 



D + q 2 B + g 4 C + q 6 E , D 1 + q 2 B 1 + q 4 C 1 



+ 



D-7 



U 2 V 3 

q 2 B 2 , D 3 



U 3 V 2 



U 4 V 



Vo 
V 



(19) 



when the D-operation is carried out (omitting the factor (a/V) 4 for simplicity). The co- 
efficients Di, Bi, Ci, and E\ are functions of building blocks, Bij and Aj described in the 
following subsection. The suffix I labels the number of contractions of D-operators. 



D. Building blocks 

Building blocks, U and B^, are homogeneous polynomials of degree n and n — 1 of the 
Feynman parameters Zi, respectively. They are determined by the topological structure of 
the diagram called the chain diagram that is derived by amputating all external lines and 



disregarding distinction of the types of lines [25] 
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The fundamental set of circuits of the chain diagram that consists of n independent 
self-nonintersecting closed loops are chosen in the following manner. One type of circuit is 
formed by an internal photon line and consecutive lepton lines that connect the endpoints 
of the photon line. The direction of the circuit is chosen to be that of the lepton lines. We 
may assign the direction of the photon line accordingly. There are n — 1 circuits of this sort. 
The n-th circuit is chosen to be the closed lepton loop itself. 

Then, for i and j that label the lines of chain diagram, U and Bij are given by 



U st = y~] z k £ k ,s£k,t, U = det U st , 

' * st 

k 

B H = V (20) 

s,t 

where s and t refer to the circuits. The loop matrix £k, c takes (1, —1, 0) according to whether 
the line i is (along, against, outside of) the circuit c. 

Once U and are obtained, another building block, the scalar current Aj, is given by 

A i = Yl (fy-ZiBy/U). (21) 

E. UV divergence 

The amplitude constructed thus far may have UV divergences when the sum of Feynman 
parameters of one or more internal loops tends to zero. We adopt subtractive renormalization 
here in a suitable way for numerical treatment. The subtraction terms are prepared as an 
integral over the same Feynman parameter space as the original unrenormalized amplitude 
so that the divergences of the amplitude are canceled pointwise. These subtraction terms are 
constructed by a simple algorithm called i^-operation [25( for each occurrence of subdiagram 
UV divergences. The whole divergent structure of a diagram is recognized by Zimmermann's 
forest formula. 

The subdiagrams relevant for the UV divergence are of self-energy type or vertex type. 
For the proper lepton loops of the present concern, such a subdiagram is represented by an 
open segment of the lepton loop, which involves vertices and lepton lines in the segment, 
and photon lines whose endpoints are included in the segment. Therefore, we have to exam- 
ine every segment of the lepton loop whenever it corresponds to a one-particle irreducible 
subdiagram of self-energy type or vertex type. 
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The inclusion relation of subdiagrams are found by examining the inclusion relation of 
the segments: they are independent or overlapping, or one segment is completely included in 
the other. Once the relation is known, the complete set of forests are constructed by finding 
all possible sets of subdiagrams whose elements are not overlapping with each other. 

Each forest corresponds to a particular emergence of UV divergence, and it is related 
to a subtraction term in the subtractive renormalization. For a vertex subdiagram S of 
a diagram Q, the subtraction term defined by the i^-operation factorizes analytically by 
construction as 

^ v n s/5 , (22) 

where L^ v is the UV-divergent part of vertex renormalization constant for the subdiagram 
S, and Hg/s is the amplitude of the reduced diagram Q/S. When a subdiagram S is of 
self-energy type, the subtraction term factorizes analytically as 

5m^U g/sin + B^U g/Stl/ , (23) 

where 5m^ v is the UV-divergent part of mass renormalization constant 8ms, and Bg is the 
UV-divergent part of wave-function renormalization constant B$- When a forest consists 
of more than one subdiagram, the subtraction term becomes products of renormalization 
constants and reduced amplitudes. 



III. INSERTION OF n< 8 ' IN M 2 



The easiest way to insert the eighth-order vacuum-polarization loops in M 2 is by the 
formula 



29 



30] 



M> 



2,P 8 



V 1 - V 



(24) 



where II® is given by Eqs. f fl8|) and f fl9|) . It is straightforward to include this in the 
automated code generation system. 

As a check of the integration codes, we have also derived the following formula from 
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Eq. (|T9|) applying the method described in Ref. |25[ which yields 
-i 



M. 



2,P& 



where 



dy(l - y) / (dz) 
Co 



D ° (w 3 -(W- if) + — 

UW*W 3K K ' ' U 2 GV 2 W 3 



U 2 G 2 V W 3 
B 1 



(W-l) + 



E 



+ 



U 2 G 3 W 3 U*V 2 W 2 



(W 2 -{W - if) 



U 3 GVoW 2 



(W-l) 



+ 



D 7 



B 2 



In 



U 3 G 2 W 2 U A V a W U A GW U 5 \W-1 



(W-l) 



w 



w 



V l-y 



(25) 



(26) 



G y 2 

The UV divergent terms of M 2j p g can be isolated by the i^-operation. See Appendix O for 
details. 



IV. RESIDUAL RENORMALIZATION 

The standard on-the-mass-shell renormalization of our vacuum-polarization function IT® 
is given explicitly in Appendix [B] Actually, it is not suitable for evaluation of these terms 
on the computer, because individual terms of these functions are UV-divergent. Thus it is 
necessary to convert them into sums of UV-finite quantities. This is achieved by an interme- 
diate renormalization procedure carried out by the i^-operation shown in Appendix O The 
i^-operation subtracts only the UV-divergent parts of renormalization constants. In order to 
obtain the standard on-the-mass-shell renormalization, the remaining UV-finite terms must 
be removed by a procedure called residual renormalization. 

Substituting expressions given in Appendix O into corresponding expressions in Ap- 
pendix El and making use of various subdiagram relations listed in Appendix [Dj we can 
convert the right-hand-side of equations of Appendix [B] into the sum of finite quantities. 
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Collecting all terms thus created we obtain 
a (lO) [I(i) (M 2 ) ] = M ( S 

- 6ALB 2 M«$ 6 

+ {14(ALB 2 ) 2 - 4ALS 4 }M 2 ( ;g 

+ {-14(ALB 2 ) 3 + 14ALS 4 ALB 2 - 2ALB e }M^ff 

- A5m A M^l 

+ (12ALB 2 A5m 4 + 2A5m 4 A5m* 2 - 2A5m 6 )M^. (27) 

Suppressing the superscript (hh) for simplicity the residual renormalizaton terms are defined 
as follows: 

p39 

M 2:A p 8 = ^2 n FiM 2ti , 
i=p01 
H 

M 2iAP4 = M 2jP4a +2M 2:Pib , 
M 2iA p 4 , = M 2i p 4o ,+2M 2i p 4i „ 

ALB 6 = XpALBep, 

3 

ALB 4 = ^(AL 4tM + AL 4M ) + AE 4a + A5 4fe , 

i=l 

ALB 2 = AB 2 

i? 

A5m6 = \pA5rriGp, 

(3=A 

A<5m 4 = A<5m 4a + A<5m 4fe , (28) 

where npj = 1 for % = p01,p02,p03, npi = 2 for i = p04, ...,p22, and npi = 4 for i = 
p23, ...,p39, r] A = 7] C = Vd = Vf = 2, Vb = Vg = Vh = 1, = 4, and A A = A B = A c = 
As = A_p = A/r = 1, Ax> = Ag = 2. ALB 6/3 , (3 = A, . . . ,H, are defined in Appendix ID1 
Numerical values of ALB 2 , ASm^, A<5m 4 , M 2 .p 2 „, 5m 2 *, etc., are listed in Table HT1 

The coefficient —6 of M 2i ap 6 in (127|) can be readily understood noting that the vacuum- 
polarization function Il 6 has 6 fermion lines into which two-point vertex can be inserted. 
This insertion is the source of the wave-function renormalization constant B 2 and the self- 
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mass 5m2 term. Since A<5m2 = 0, however, only the AB2 survives in ( |27|) . For convenience 
let us call this an insertion of B2. 

We find 14 different ways of insertion of two B 2 's in II4. Ten comes from insertions of 
two disconnected second-order self-energy diagrams and four comes from insertions of the 
fourth-order self-energy diagram in which a second-order self-energy diagram is completely 
included in another second-order self-energy diagram. Three B2S can be inserted in II2 in 
14 ways. Insertion of one B4 and one B 2 in II2 can also be made in 14 different ways. All 
these three terms should be accompanied by terms proportional to A5m 2 which, however, 
vanish in our formulation based on the .fT-operation. 

The factor —4 in — 4ALB 4 M 2! ap 4 is due to the fact that II4 has 4 fermion lines into which 
-B4 can be inserted. The apparent absence of the coefficient 4 in —A5m 4: M2 t AP 4 * can be 
accounted for by the fact that II4 has four fermion lines into which a two-point vertex can 
be inserted. Thus the coefficient 4 is absorbed in the definition of n 4 » . This term is present 
in f )27|) since A<5m 4 is nonvanishing. 

Finally one B$ can be inserted in II2 in two ways. The factor 2 in —2A5rriQM2 y p 2t is 
the same as that of — 2AB§M2,p 2 . The term 2A5rri4 : A8m2*M2,p 2t is related to the sub- 
diagram of — 2A<5m 6 M 2) p 2 , except for the factor —1. Application of i^-operation on the 
second-order self-energy subdiagram of 1I4* yields 6A52II2*. Application of i^-operation 
on the second-order self-energy subdiagram of 5rriQ yields 6 A5777.4A.B2. Together they give 
12ALB2A5m 4 M 2i p 2< ,. 

A similar argument can be given starting from vertex renormalization subdiagrams al- 
though it does not give information on A5m term. Consideration on vertex renormalization 
is not necessary, however, since a[ [I(i)] is free from infrared(IR) divergence so that L n is 
always combined with B n to form an finite combination ALB n . 

The reason the coefficients of residual renormalization terms can be determined by the 
argument described above is that the UV-finite parts are not affected by i^-operation which 
deals only with UV-divergent parts so that the coefficients of residual renormalization terms 
inherit the structure of the standard renormalization unaltered. 
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TABLE I: Contributions of diagrams of Set I(i) to a e for (l\l 2 ) = (ee). The superscript (ee) is 
suppressed for simplicity. The multiplicity tif is the number of vertex diagrams represented by the 
integral and is incorporated in the numerical value. All integrals are evaluated initially with 10 8 
sampling points per iteration, iterated 150 times, followed by 10 9 points, iterated 10 times. 



Integral 


np 


Value (Error^ 
including rip 






Sampling per 
iteration 


No. of 
iterations 


M 2 , p 01 


1 


0.035 


760 8 


(60) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 




1 


0.017 


303 9 


(40) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2 , P 03 


1 


0.039 


757 3 


(73) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


-M"2,p04 


2 


0.037 


755 2 


(93) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


-^2,p05 


2 


0.062 


796 


(190) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


A^2,p06 


2 


0.129 


748 


(225) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2 , P 07 


2 


0.128 


655 4 


(251) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M2, P 08 


2 


0.103 


304 6 


(167) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M2,p09 


2 


-0.038 


968 7 


(92) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


A^2,pl0 


2 


-0.057 


281 7 


(105) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2 , P n 


2 


0.020 


893 9 


(34) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2 ,pl2 


2 


0.038 


800 8 


(54) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2 , P 13 


2 


0.017 


581 6 


(166) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2 ,pl4 


2 


0.090 


813 8 


(165) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 




2 


0.008 


522 3 


(17) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2iP i 6 


2 


0.023 


211 8 


(28) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2tPl7 


2 


-0.009 


190 1 


(34) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2 , p i8 


2 


0.011 


705 8 


(31) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2 ,pi 9 


2 


0.024 


548 


(48) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2jP 20 


2 


0.028 


129 1 


(24) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2jp2 i 


2 


-0.014 


422 7 


(29) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2tp22 


2 


0.012 


984 7 


(39) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2 , p23 


4 


0.008 


650 4 


(35) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2)P2i 


4 


0.037 


873 8 


(95) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2jP 25 


4 


0.168 


761 3 


(302) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2jp26 


4 


-0.061 


567 5 


(141) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2tP27 


4 


-0.145 


657 6 


(325) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2)P2S 


4 


-0.078 


819 8 


(148) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2jp29 


4 


0.110 


765 8 


(297) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2 ,p 30 


4 


0.217 


591 7 


(407) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2jp31 


4 


-0.149 


396 6 


(349) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2tP32 


4 


-0.122 


439 


(139) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2yp33 


4 


0.043 


600 9 


(107) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2ip3 4 


4 


-0.003 


177 4 


(75) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2 , p35 


4 


-0.054 


641 9 


(129) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2jp3 e 


4 


-0.138 


680 8 


(150) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2 , p37 


4 


-0.102 


260 8 


(164) 


1 


X 


10 s , 


1 


X 


10 9 


150, 10 


M 2)P3S 


4 


-0.060 


892 7 


(234) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 


M 2:P39 


4 


-0.008 


689 5 


(43) 


1 


X 


10 8 , 


1 


X 


10 9 


150, 10 
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TABLE II: Auxiliary integrals for Set I(i). Some integrals are known exactly. Other integrals are 
obtained by the integration routine VEGAS. The suserscript (li I2) indicates that the open and 
closed fermion lines consist of fermions l\ and I2, respectively. The letters e, m, and t stand for 
electron, muon, and tau-lepton, respectively. 



T 4- 1 

Integral 


Value (error) 


Integral 


Value (error) 


A5m* 2 


-0.75 




ALB 2 


0.75 




ASm^ 


1.906 


340 (21) 


AZ£?4 


0.027 


930 (28) 


AdrriQ 


-2.340 


65 (48) 


ALBq 


0.100 


86 (77) 


1V1 2,P2 


0.015 


687 421 •• • 


M^l 


0.519 


762 (21) x 10~ 6 


M^l 


1.094 


258 28 • • • 


M^l 


0.000 


078 067 (4) 


M&l 


-0.012 


702 383 •• • 


mM 

1V1 2,P2* 


-0.519 


719 (17) x 10~ 6 


M (me) 


-0.161 


084 05 • • ■ 


M (mt) 
IV1 2,P2* 


-0.000 


077 655 (3) 


M (ee) 

JW 2,AP4 


0.076 


401 785 •• • 


mM 

iw 2,AP4 


0.275 


271 (1) x 10~ 5 


M (me) 
JW 2,AP4 


3.135 


059 01 (2) 


M (mt) 
JW 2,AP4 


0.000 


412 61 (3) 


M (ee) 

JW 2,AP4* 


-0.117 


770 (12) 


M (em) 
iw 2,AP4* 


-0.000 


005 505 (1) 


M (me) 
iw 2,AP4* 


-0.754 


40 (13) 


M (mt) 
iw 2,AP4* 


-0.000 


819 49 (8) 


iw 2,AP6 


0.187 


046 (7) 


M {em) 
JW 2,AP6 


0.731 


632 (71) x 10~ 5 


M (me) 
iw 2,AP6 


5.543 


94 (42) 


M (mt) 
iw 2,AP6 


0.001 


094 28 (7) 



V. NUMERICAL EVALUATION OF M 



(ee) 



FORTRAN codes of our diagrams are generated by gencodevpA following the proce- 
dures described in Section [III The validity of gencodevp A has been tested thoroughly for 
diagrams of Set 11(d) whose integrals are known by several other means Numeri- 
cal integration was carried out by an adaptive-iterative Monte-Carlo numerical integration 



routine VEGAS 



35] with 10 8 sampling points per iteration and 150 iterations followed by 



10 9 sampling points per iteration and 10 iterations. The results are summarized in Tables [I] 
and [III From these data we obtain 



4 10) [i(i 



,(ee)l 



0.017 47 (11). 



(29) 



VI. NUMERICAL EVALUATION OF 

Once FORTRAN programs for mass-independent diagrams are obtained, it is straightfor- 
ward to evaluate the contribution of mass-dependent term A^^e/m^). We simply have 
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TABLE III: Contributions of diagrams of Set I(i) to a e for (Z1/2) = (em). The multiplicity np is the 
number of vertex diagrams represented by the integral and is incorporated in the numerical value. 
The superscript (em) is omitted for simplicity. All integrals are evaluated in double precision. 



Integral 


np 


Value (Error) 
including np 


Sampling per 
No. of iteration 


No. of 
iterations 


M 2 , p oi 


1 


0.115 63 


(3) x 10~ 5 


1 X 


10 8 


50 


M 2 , p0 2 


1 


0.077 01 


(3) x 10~ 5 


1 X 


10 8 


50 


M 2 , P 03 


1 


0.201 24 


(5) x 10~ 5 


1 X 


10 8 


50 


-^2,p04 


2 


0.147 75 


(6) x 10~ 5 


1 X 


10 8 


50 


M2,p05 


2 


0.236 81 


(14) x 10~ 5 


1 X 


10 8 


50 


M2,p06 


2 


0.484 29 


(16) x 10~ 5 


1 X 


10 8 


50 


M2,p07 


2 


0.488 27 


(19) x 10~ 5 


1 X 


10 8 


50 


M2, P 08 


2 


0.400 43 


(12) x 10~ 5 


1 X 


10 8 


50 


M2, P 09 


2 


-0.148 04 


(6) x 10~ 5 


1 X 


10 8 


50 


M 2 , P 10 


2 


-0.248 34 


(8) x 10~ 5 


1 X 


10 8 


50 


M 2 , P n 


2 


0.071 31 


(2) x 10~ 5 


1 X 


10 8 


50 


M 2 ,pl2 


2 


0.168 94 


(4) x 10~ 5 


1 X 


10 8 


50 


M 2 , P 13 


2 


0.052 85 


(12) x 10~ 5 


1 X 


10 8 


50 


M 2 ,pl4 


2 


0.409 41 


(12) x 10- 5 


1 X 


10 8 


50 




2 


0.043 04 


(1) x 10- 5 


1 X 


10 8 


50 


M 2 , P w 


2 


0.064 53 


(1) x 10" 5 


1 X 


10 8 


50 


M 2 , P n 


2 


-0.044 76 


(2) x 10" 5 


1 X 


10 8 


50 


A^2,pl8 


2 


0.037 80 


(2) x 10" 5 


1 X 


10 8 


50 


M 2 , P 19 


2 


0.088 50 


(3) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p20 


2 


0.114 22 


(1) x 10" 5 


1 X 


10 8 


50 


M 2 ,p21 


2 


-0.062 59 


(2) x 10~ 5 


1 X 


10 8 


50 


M2,p22 


2 


0.047 81 


(2) x 10~ 5 


1 X 


10 8 


50 


M 2 , P 23 


4 


0.034 72 


(2) x 10" 5 


1 X 


10 8 


50 


M 2 , P 24 


4 


0.148 88 


(7) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p25 


4 


0.652 31 


(23) x 10" 5 


1 X 


10 8 


50 


M 2 ,p26 


4 


-0.235 83 


(10) x 10" 5 


1 X 


10 8 


50 


M2,p27 


4 


-0.491 34 


(21) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p28 


4 


-0.347 15 


(11) x 10" 5 


1 X 


10 8 


50 


M 2 ,p29 


4 


0.389 88 


(21) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p30 


4 


0.944 16 


(33) x 10- 5 


1 X 


10 8 


50 


M 2 , P 31 


4 


-0.512 63 


(25) x 10~ 5 


1 X 


10 8 


50 


M2,p32 


4 


-0.477 71 


(9) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p33 


4 


0.146 82 


(6) x 10" 5 


1 X 


10 8 


50 


M 2 , P 34 


4 


-0.000 10 


(5) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p35 


4 


-0.171 64 


(9) x 10- 5 


1 X 


10 8 


50 


M 2 ,p36 


4 


-0.548 36 


(11) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p37 


4 


-0.498 39 


(13) x 10~ 5 


1 X 


10 8 


50 


M 2 ,p38 


4 


-0.179 07 


(15) x 10" 5 


1 X 


10 8 


50 


M 2 ,p39 


4 


-0.044 15 


(3) x 10~ 5 


1 X 


10 8 


50 
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TABLE IV: Contributions of diagrams of Set I(i) to for (l\l 2 ) = (me). The multiplicity np is 
the number of vertex diagrams represented by the integral and is incorporated in the numerical 
value. The superscript (me) is omitted for simplicity. All integrals are evaluated initially with 
sampling points 10 8 per iteration, iterated 50 times, followed by 10 9 points per iteration, iterated 
200 times, and 10 10 points, iterated 10 to 80 times. 



Integral np Value (Error) Sampling per No. of 

including np iteration iterations 



M 2 , P oi 


1 


5.475 765 


(50) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


40 


M 2 , P 02 


1 


-3.639 035 


(13) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


40 


M 2 , P 03 


1 


1.014 582 


(56) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


40 




2 


9.957 281 


(77) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M2, P 05 


2 


11.130 636 


(109) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


80 


M2, P 06 


2 


3.445 706 


(95) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M 2 , P 07 


2 


1.150 328 


(136) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


70 


M2,p08 


2 


2.431 621 


(98) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M2, P 09 


2 


-6.305 904 


(87) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M2, P 10 


2 


3.576 267 


(69) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


10 


M 2 , p u 


2 


3.087 991 


(41) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 ,pl2 


2 


2.681 026 


(32) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


10 


M 2 , P 13 


2 


8.166 698 


(70) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 ,pl4 


2 


-2.042 862 


(97) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M 2 ,pl 5 


2 


-0.014 213 


(9) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 , p w 


2 


3.556 341 


(11) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 , P n 


2 


3.279 641 


(8) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 , P 18 


2 


-1.044 365 


(5) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M2, P 19 


2 


3.238 053 


(31) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M2,p20 


2 


-0.804 464 


(15) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 , P 21 


2 


-2.002 941 


(28) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 ,p22 


2 


1.439 287 


(14) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 ,p23 


4 


2.246 387 


(9) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 




4 


-1.943 020 


(59) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


50 


M 2 ,p25 


4 


0.717 756 


(211) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


80 


M 2 ,p26 


4 


-11.597 221 


(131) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


80 


M2,p27 


/I 


1 R AQ7 QQfl 

— 10. 4y/ yyu 




1 
1 


X 


in 8 

1U , 


1 

1 


X 


in 9 

1U , 


1 
1 


X 


1U 


ou, 


zuu, 


an 

oU 


M 2 , P 28 


4 


7.463 932 


(85) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M 2 , P 29 


4 


18.659 291 


(133) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


80 


M2,p30 


4 


-3.240 940 


(198) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M 2 , P 31 


4 


-16.369 751 


(181) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


80 


M 2 , P 32 


4 


-1.501 462 


(99) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


70 


M 2 , P 33 


4 


8.005 147 


(70) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M 2 , P 34 


4 


-5.994 355 


(74) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M2,p35 


4 


-13.319 800 


(61) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 


M2,p36 


4 


6.372 372 


(64) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


50 


M 2 ,p37 


4 


-0.889 232 


(94) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M 2 ,p38 


4 


-12.751 250 


(97) 


1 


X 


10 s , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


60 


M 2 ,p39 


4 


1.541 038 


(33) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


200, 


20 
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TABLE V: Contributions of diagrams of Set l(i) to for (hfo) = (mi). The multiplicity rip is the 
number of vertex diagrams represented by the integral and is incorporated in the numerical value. 
The superscript (mi) is omitted for simplicity. 



Integral 


np 


Value (Error) 
including 




Sampling per 
No. of iteration 


No. of 
iterations 




1 


0.000 173 73 


(5) 


1 x 10 8 


50 


M 2)p0 2 


1 


0.000 115 02 


(4) 


1 x 10 8 


50 


M 2)P 03 


1 


0.000 299 06 


(8) 


1 x 10 8 


50 


M 2 , p 04 


2 


0.000 220 89 


(9) 


1 x 10 8 


50 


M2,p05 


2 


0.000 354 37 


(21) 


1 x 10 8 


50 


M2,p06 


2 


0.000 725 06 


(24) 


1 x 10 8 


50 


M2,p07 


2 


0.000 730 62 


(28) 


1 x 10 8 


50 


M2,p08 


2 


0.000 599 02 


(19) 


1 x 10 8 


50 


M2,p09 


2 


-0.000 221 58 


(10) 


1 x 10 8 


50 


M 2 ,p W 


2 


-0.000 370 54 


(12) 


1 x 10 8 


50 


M 2 , p u 


2 


0.000 107 04 


(3) 


1 x 10 8 


50 


M 2 ,pl2 


2 


0.000 252 11 


(6) 


1 x 10 8 


50 


M 2 ,pl3 


2 


0.000 079 35 


(18) 


1 x 10 8 


50 


M 2 ,plA 


2 


0.000 610 34 


(18) 


1 x 10 8 


50 


-^2,pl5 


2 


0.000 064 07 


(1) 


1 x 10 s 


50 


M 2 ,pW 


2 


0.000 097 35 


(2) 


1 x 10 s 


50 


M 2 , P n 


2 


-0.000 066 71 


(3) 


1 x 10 8 


50 


M 2 ,pis 


2 


0.000 056 95 


(3) 


1 x 10 8 


50 


M 2 , P 19 


2 


0.000 132 63 


(4) 


1 x 10 8 


50 


M 2 ,p2Q 


2 


0.000 170 72 


(2) 


1 x 10 8 


50 


M 2 ,p21 


2 


-0.000 093 37 


(3) 


1 x 10 8 


50 


M2 tP 22 


2 


0.000 071 59 


(4) 


1 x 10 8 


50 


M 2 ,p23 


4 


0.000 051 80 


(4) 


1 x 10 8 


50 


M 2 ,p2A 


4 


0.000 222 52 


(11) 


1 x 10 8 


50 


M 2 ,p25 


4 


0.000 975 69 


(35) 


1 x 10 8 


50 


M 2 ,p26 


4 


-0.000 352 91 


(15) 


1 x 10 8 


50 


M2, P 27 


4 


-0.000 737 37 


(32) 


1 x 10 8 


50 


M 2 ,p2S 


4 


-0.000 517 81 


(17) 


1 x 10 8 


50 


M 2 ,p29 


4 


0.000 584 31 


(32) 


1 x 10 8 


50 


M 2 ,p30 


4 


0.001 408 74 


(49) 


1 x 10 8 


50 


M 2 , P 31 


4 


-0.000 768 96 


(38) 


1 x 10 8 


50 


M 2 , P 32 


4 


-0.000 714 49 


(14) 


1 x 10 8 


50 


M 2 , P 33 


4 


0.000 220 34 


(10) 


1 x 10 8 


50 


M 2 , P 3A 


4 


-0.000 000 59 


(8) 


1 x 10 8 


50 


M 2 ,p35 


4 


-0.000 257 99 


(13) 


1 x 10 s 


50 


M 2 ,p36 


4 


-0.000 820 21 


(17) 


1 x 10 s 


50 


M 2 ,p37 


4 


-0.000 741 52 


(19) 


1 x 10 s 


50 


M 2 ,p38 


4 


-0.000 269 58 


(22) 


1 x 10 s 


50 


M 2 ,p39 


4 


-0.000 065 68 


(4) 


1 x 10 8 


50 
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to choose an appropriate loop fermion mass. The results are summarized in Table IIHI From 
this table we obtain 

a( 10) [I(i) W ] = 0.000 001 666 (24). (30) 

The numerical data used to obtain ( 1301) are listed in Tables IIHI and [Til The contribution 
A^2°\m e /m T ) is two orders of magnitude smaller than (1301) so that it is negligible at present. 



VII. CONTRIBUTION TO THE MUON g-2 

The codes described above can also be applied to calculate the contribution of the Set I(i) 
to the muon g — 2. From Tables ITVl and HTl we obtain 



( 10 )[I(i) (me) ] = 0.087 1 (59). (31) 



We have also evaluated the tau-lepton contribution. The values listed in Tables [V] and [TTJ 
lead to 

a^[l(i) (mt) ] = 0.000 237 (1). (32) 
The contribution of Set I(i) diagrams to muon g — 2 was first discussed in Eq. (19) of 



Ref. 



26] in which the renormalization group was efficiently used to pick up the leading 



logarithmic contribution: 



oni(i)M] = _i afl + m _ a, ln ^ + Q(mJm ^ (33) 

where is the constant term of the asymptotic expansion of the proper vacuum- 
polarization function II® (q) of the eighth order, which was left undetermined. This constant 
can be determined from our numerical result Eq. ( I3TI) : 



af 1 [numerical] = -1.641 (59) . (34) 

Recently the asymptotic analytic form of afi was obtained directly together with other 
eighth-order vacuum-polarization diagrams [281] . The explicit form of their a^ 1 ' was given 



only in the slide of the conference talk [27J, which reads 

ajj [anal. -asympt.] = —1.971 6 +0(m e /m^). (35) 
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Substituting this value to Eq. ( 133|) . they obtained the asymptotic contribution of the Set I(i) 

aj t 10) [I(i) : anal.-asympt.] = 0.252 37 + 0(rn e /m M ). (36) 

Whether this is in agreement with our result (1311) or not is somewhat subtle and will be 
discussed in the next section. 

VIII. SUMMARY 

In this paper we obtained the eighth-order vacuum-polarization function IT( 8 )(g 2 ) as a 
sum of Feynman-parametric integrals. It is then applied to the calculation of the tenth-order 
lepton g — 2. Collecting (129]) and fl30l) we obtain the contribution of the gauge-invariant Set 
I(i) to the electron g — 2 

c4 10) [I(i) (all) ] = 0.017 47 (11). (37) 
From (I2"9l) . (|3"Tj) and (1321) we obtain the contribution from Set I(i) to the muon g — 2 

a^[I(i) (all) ] =0.104 8 (59). (38) 

It is difficult to decide whether our result (|3T|) and asymptotic result (|36|) are in agreement 
or not. In order to illuminate this problem it may be helpful to compare of similar 
structure in lower orders. 

For the 6th-order a M obtained by inserting a proper 4th-order vacuum-polarization IT^ 



in the second order M2 gives 
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36| 



aJfVum.] = 1.493 671 581 (8), 

ajf^asym.] = 1.517 3 + 0(m e /m M ), (39) 

where the overall factor (a/n) 3 is omitted for simplicity. 

For the 8th-order a M obtained by inserting a proper 6th-order vacuum-polarization II^ 
in the second order M 2 gives the coefficients of (a/V) 4 : 

[num.] = -0.230 596 (416), 
ajfpade] = -0.230 362 (5), 

ajf^asym.] = -0.290 987 + 0(m e /m M ), (40) 
22 



is from 



[EL 



where the numerical evaluation aj?[num.] is from [331, |37| . the Pade approximation a® [Pade] 



3l|, and the asymptotic result a®[asym.] is from [38|. Note that the asymptotic 
result of fHUj) contains the leading logarithmic and next-to-leading constant terms. 

The difference between the numerical and asymptotic results come from the contribution 
of order m e /m^. From the sixth-order ( )39|) . eighth-order ( 140]) . and tenth-order cases we find 

af [num.] - ajf } [asym.] = -0.024, 
af [num.] - af [asym.] = 0.061, 

a£°Vum.] - ^ 10) [ as y m -] = °- 17 • ( 41 ) 

Note that the difference increases as the order of perturbation increases. Nevertheless, we 
cannot exclude the possibility that the difference between ( 13~TT) and ( 1361) is caused by some 
error. One possible cause is that the uncertainty of ( )3T|) is gross underestimate because of 
insufficient data sampling. The situation might be similar to the case of aj, [num.] in early 



calculations 
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39 



401 ] where poor sampling of integrands was found to be the cause of 
large discrepancy with the Pade result. This problem was finally settled by going to a much 
larger sampling, which led to (|40|) . 

In order to examine the possibility of gross underestimate of errors in ( l3Tj) we evaluated 
the integrals with the sampling points per iteration M of 10 8 , 10 9 , and even with 10 10 . 
The results show no sign of the central values drifting beyond the error bars estimated by 
VEGAS as M increases. We are therefore confident that our result fl3T|) is free from the 
problem caused by insufficient samplings. 
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(a)n 4a (b)n 4b 

FIG. 5: Fourth-order vacuum polarization diagrams. 




(a)n 4a ,i» (b)n 4M . (c)n 4M * (d)n 4fei4 * 

FIG. 6: Fourth-order vacuum polarization diagrams with mass insertions. 

Appendix A: Fourth-order vacuum-polarization functions with mass insertion and 
their contribution to M2 p 4 * 

Since M 2) ap 4 » does not appear except in Set I(i) in our study of the tenth-order g—2, we 
shall derive parametric formulas for the diagrams n 4a * , U ib * , and M 2t p it in this Appendix. 

Our derivation follows closely the treatment of diagrams and IT^ (without mass 
insertion) which consist of four lepton lines forming a closed loop and an internal photon 
line a as shown in Fig. [5j Following the steps leading to Eq. ( fl4"j) of Sec. Ill CI we obtain 




From Lorentz invariance and gauge invariance, we have the general structure 

= {q^qu - q^gtiu)^* + (gauge-dependent terms). (A2) 
Charge renormalization is achieved by 

n 4 ,(g 2 ) = n 4 »(g 2 ) -n 4 .(o). (A3) 
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When the D-operation is carried out in il 4 * , the result can be expressed in the form 



(dz) 



Do 



U 2 [ V 2 V,? J + V 2 



q 2 B , D 



f U 3 \V Vq 



where Vq = Zi2um 2 and D , B , D\, U, V and (dz) are diagram-specific. 



(A4) 



1. Diagram M 2> p Aait 

This diagram has four fermion lines into which mass vertex can be inserted. They all 
give the same contribution to M 2) p 4a » so that we have to evaluate only one of them such as 



n 4a l » of Fig. 6(a) For this diagram we find 



D = r\-2A 3 A A - 2A 2 A 4 + A 2 A 3 + 2A X A A - A X A 3 - A X A 2 ), 

Bo = r 2 A\{A 2 A 3 -A l A 3 -A x A 2 ), 

D l = r 2 B 11 (A 3 A 4 + A 2 A i + 2A 1 A 3 + 2A 1 A 2 ) 

+ r 2 B 12 (-3A 3 A A - A 2 A A + 2A t A 4 - 2A l A 2 ) 

+ r 2 B 22 (2A 1 A 4 ) 



(AS) 



where r is the ratio of the mass m of the loop lepton and the mass of the lepton of M 2 and 



^3 

B n 
U 
(dz) 



1 - ZiBn - z 2 B 12 , A 2 = l- ziB 12 - z 2 B 22 , 



-z x B l3 - z 2 B- 



23; 



A A 



-z-lB u - z 2 B- 



24- 



^23a; B\2 — Z a , B 22 — Zua, 

z Ua Z23 + z u z a , G = z 1 A 1 + z 2 A 2 , V = V - q 2 G, 
z 1 dz 1 dz 2 dz 3 dz A dz a 8(l - z 123Aa ), z % > 0. 



(A6) 



This diagram has a UV divergence from the subvertex {2, 3, a,}, which can be isolated by 
the i^23 operation. Subtraction of this term yields a UV-finite value 



An 4a * = 4(1 - K 23 )U 4ai 



1*, 



(A7) 



and a finite contribution M 2y p iat , to g — 2. By numerical integration we obtain the value 
-0.066 907 (7) for the (ee) case, and -0.443 935 (69) for the (me) case. 
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2. Diagram M 2 p.., 



The diagrams IT^i* and n 4 b j3 * give identical contribution to M 2t p 4b , , whereas the diagrams 
n 4fe2 * and n 4b4 * must be treated separately. 



For the diagrams n 461 . (see Fig. 6(b) ) we find 



D = r\2A 2 A 4 -QA 1 A A ), 
B = r 2 (3A 2 1 A 2 A 4 -2A 3 1 A 4 ), 
D x = r 2 B 11 (6A 1 A 4 -A 1 A 2 + 2A 2 l ) 
+ r 2 B 12 (-9A 1 A 4 -2A 2 1 ), 



(A8) 



where 



At 
A 2 
A 3 
B u 
B\ 2 
U 
(dz) 



1 - z x B X \ - z 2 B l2 - z 3 B 13 , 

1 — z\B\ 2 — z 2 B 22 — z 3 B 23 , 

A x , A A = A 1 - 1, 

B13 = B14 = B 33 = B 3 4 = B 44: = z 2a , 

B 23 = B 24 = Za, B 22 = Z\ 3 ia, 

ziuz 2a + z 2 z a , G = z x A x + z 2 A 2 + z 3 A 3 , V = V - q 2 G, 
zidzidz 2 dz 3 dz 4 dz a 5(l - z l234a ), z { > 0. 



(A9) 



This diagram has a UV-divergence from the self-energy subdiagram {2, a} which can be 
subtracted by the ^-operation 



AIT 4bil * = (1 - K 2 )U 



4b,l* 



(A10) 



For the diagrams n 46 2 * (see Fig. 6(c)) we find 



Do = r\2A 2 A 4 -AA 1 A A ), 
B = r 2 (-AA l A 2 2 A 4 + 2A 2 1 A 2 A 4 ), 
Dx = r 2 B 22 (8AxA 4 ) 
+ r 2 B 11 {2A 2 A 4 ) 

+ r 2 B 12 (AA 2 A 4 - 8AxA 4 + 4AxA 2 - 2A\), 



(All) 
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where 

At = 1 - Zi-Bn - 2 2J Bi2 - 2 3j B 13 , A 2 = 1 - 2iSia - z 2j B 22 - z 3j B 23 , 

A 3 = 1 - Zi-B 13 - z 2 B 23 - z 3 B 33 , 

A4 = —ziBu — z 2 B 2 4 — Z3B34, 

Bn — B iS = B14 = B 33 = B 34 = B44 = z 2ai 

Byi = B 23 = B 24 = z a , B 22 = 2 134a , 

U = z lu z 2a + z 2 z a , G = z x A x + z 2 A 2 + z 3 A 3 , V = V - q 2 G, 

(dz) = z 2 dzidz 2 dz 3 dzidz a 5{l — Zx 23 4 a ), Zj > 0. (A12) 

This diagram has a UV- divergence from the subdiagram {2*, a} which can be subtracted 
by the K 2 * -operation 

An 46 , 2 * = (1-K 2 *) n 46,2*. (A13) 



For the diagrams n 4bi4 * (see Fig. 6(d)) we find 

D = r\2A 2 A 4 -8A 1 A 4 ), 
B = r 2 (2AjA 2 A 4 ), 
D 1 = r 2 Bxi{2A 2 A A ) 

+ r 2 B 12 (-8A 1 A 4 -2A 2 1 ), (A14) 

where 

Ax = 1 - ZxB n - z 2 B 12 - z 3 B 13 , A 2 = 1 - z x B 12 - z 2 B 22 - z 3 B 23 , 
A 3 = 1 - zxB 13 - z 2 B 23 - z 3 B 33 , 
A 4 = —ZxBx4 — z 2 B 24 — z 3 B 34 , 
Bxx = Bx3 = Bxa = B 33 = B 34 = B 44 = z 2a , 

Bx2 — B 23 = B 24 = Z a , B 22 = Z l3 4cn 

U = z 134 z 2a + z 2 z a , G = zxAx + z 2 A 2 + z 3 A 3 , V = V - q 2 G, 
(dz) = z 4 dzxdz 2 dz 3 dz 4 dz a 5{l — zx 234a ), %i > 0. (A15) 

As is for the diagram II^i* this diagram has a UV-divergence from the self-energy sub- 
diagram {2, a} which can be subtracted by the ^-operation. 
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By numerical integration we obtain for the (ee) case 

2M 2,AP ib * = -°- 046 309 ( 7 ) + °- 033 507 ( 4 ) - °- 038 061 ( 6 ) 

= -0.050 863 (9), (A16) 

where the right-hand-side of the first line is listed in order of P^i*, Pib,2*, and -P^*. 
The result for the (me) case is 

2M 2AP ib * = -°- 409 550 ( 61 ) + °- 528 759 ( 51 ) - °- 429 683 ( 62 ) 

= -0.310 464 (101). (A17) 

The sums M 2i ap 4 , = M 2 ,AP 4a , + 2M 2i AP 4b » for the (ee), (em), (me) and (mt) cases are 
listed in Table [TT1 



Appendix B: Standard on-the-mass-shell renormalization 

This Appendix describes the standard on-the-mass-shell renormalization of vacuum- 
polarization function U^ n \ where n = 2,4,6,8. n^ 2 ^ consists of only one diagram, but 
higher order functions consist of several diagrams, which must be distinguished by an ad- 
ditional symbol. For instance ]T 4i ) with i = a,b, H^' with j = A, B, H. However, the 
eighth-order functions are denoted as U^ k \ k = p01,p02, ...,p39 to avoid overcrowding. 

Renormalization terms include functions such as IT 2 *) which means insertion of a two- 
point vertex (such as a mass vertex) in the fermion line. 

Quantities L n ,B n , 5m n denote vertex renormalization constant, wave function renormal- 
ization constant, mass renormalization constant of n-th order of the standard on-the-mass- 
shell renormalization, respectively. We must also deal with renormalization constants with 
mass insertion. For instance L 2 , which contains two electron lines, it is necessary to distin- 
guish the lines into which two-point vertex insertion is made. Suppose we name them line 1 
and line 2. Then L 2 (i*i*) implies that two two-point vertices are inserted in the fermion 
line 1 of L 2 , while L 2 (i* 2 *) means that one two-point vertex is inserted in line 1 while an- 



other is inserted in line 2. (Previously 25j we used the notations L 2 **t and L 2 *t* for these 
quantities.) 
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1. Standard renormalization of fourth-order vacuum-polarization 

= n( 4fl ) - 2L 2 n< 2 > 
uifj = u^ -5m 2 u^ -b 2 u (2 ^ 

2. Standard renormalization of sixth-order vacuum-polarization 

= + 25m 2 B 2 Yi^ - 25m 2 U^ + (5m 2 ) 2 U^ - 2B 2 U^ + (B 2 ) 2 U^ 
= n< 6B > + 25m 2 B 2 Tl^ - 25m 2 U^ + (5m 2 ) 2 Ii^ - 2B 2 Tl^ + (B 2 ) 2 Ii^ 

YlfC) = n (6C) _ 5m4fc n( 2 *) + 5m 2 Jm 2 U^ + 5m 2 B 2 ,U i2) + 5m 2J B 2 n (2 * } 

- 5m 2 Ii^ - B 4b U^ - B 2 Ii^ + (B 2 ) 2 Yi^ 

n£f = n< 6D > - 5m 4a U^ + 25m 2 L 2 Tl^ - B 4a Il^ + 2B 2 L 2 U^ - 2L 2 U^ 

Ilg? = n (6E) + 6m 2 L 2 *IlW + 5m 2 L 2 U^ - 5m 2 U^ + 25 2 L 2 n( 2 ) 

- b 2 u^ - L 4M n( 2 ) - L 2 U^ 

HgP = - 2L Aa ^ - L 2 U^ + 2(L 2 ) 2 n( 2 ) 
n£f = n( 6G ) - 2L 4b , 2 n (2) - 2L 2 u^ + 3(l 2 ) 2 ii^ 
ng? = - 2L 4a , 2 n( 2 ) 
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3. Standard renormalization of eighth-order vacuum-polarization 



n (poi) = n(p oi) _ 2L 6Fi3 n( 2 ) + 4L 4a ,iL 2 n( 2 ) - 2L 2 n( 6F ) + (L 2 ) 2 U^ - 2(L 2 ) 3 n( 2 ) 

n (p02) = n ( p02 ) _ 2LeH 3U (2) 

n (po3) = n(P o3) _ 2L 6Bi3 n( 2 ) + 6L 4M L 2 n( 2 ) - 2L 4b , 2 n( 4a ) - 2L 2 n( 6G ) 
+ 3(L 2 ) 2 n( 4fl ) - 4(L 2 ) 3 n( 2 ) 

n (p04) = n ( P 04) + 65 m 2 B 2 U^ - 35m 2 (B 2 ) 2 U^ - 5m 2 U^ A ^ - 25m 2 U^ A ^ 

- 3(<5m 2 ) 2 5 2 n (2 ** ) + (<5m 2 ) 2 n (4b ' 1 * 3 * ) + 2(5m 2 ) 2 II (4M * 1 * ) - (5to 2 ) 3 II( 2 ***) - 35 2 II (6A) 

+ 3(£ 2 ) 2 n( 4fe )-( J B 2 ) 3 n( 2 ) 

n^ 5 ) = n( p05 ) + 45m 2 B 2 rt 4b ^ + 25m 2J B 2 n( 46 ' 1 *) - 35m 2 (5 2 ) 2 n( 2 *) - 25m 2 U^ B ^ 

- Sm 2 U^ A '^ - 2(5m 2 ) 2 B 2 n (2 ** ) - (5m 2 ) 2 fi 2 n (2 ** ) + (5m 2 ) 2 n (4M * 4 * } + 2(5m 2 ) 2 II (4M * 4 * ) 

- (5m 2 ) 3 n( 2 ***) - 2B 2 u^ - b 2 il<m + s(b 2 ) 2 u^ - (£ 2 ) 3 n( 2 ) 

n (p06) = n ( P 06) _ 2 5m 2 B 2 L 2 Jl^ - 25m 2 B 2 L 2 U (2 ^ + 25m 2 B 2 U^ a ^ + 25m 2 L 4b>1 U^ 

+ 25m 2 L 2 ,U ( ^ - 25m 2 U (6E ^ - 2(5m 2 ) 2 L 2 ,U (2 ^ + (5m 2 ) 2 Yl {Aa > U2 ^ + 2B 2 L 4b>1 U^ 
+ 2B 2 L 2 U^ - 2B 2 Tl^ - 2(B 2 ) 2 L 2 U^ + (B 2 ) 2 Tl^ - 2L 4bA U^ 

n ( P 07) = n ( P 07) _ 2 5m 2 B 2 L 2 Jl^ - 25m 2 B 2 L 2 U (2 ^ + 25m 2 B 2 U^ a ^ + 25m 2 L 4b>1{u) U^ 

+ 25m 2 L 2 U^ - 25m 2 rt 6E ^ - (5m 2 ) 2 L 2(U2 ^ - (5m 2 ) 2 L 2 U^ + (5m 2 ) 2 U^ U4 ^ 

+ 25 2 L 4M n (2) + 25 2 L 2 n (4b) - 2B 2 n (6B) - 2(5 2 )%n (2) + (5 2 ) 2 n (4a) 

- L 6A3 n^ - L 2 Yi^ 

n^f) = n( p08 ) - 25m 2J B 2 L 2 ,n( 2 ) - 25m 2 B 2 L 2 U {2 ^ + 25m 2 B 2 rt 4a >^ + 25m 2 L ibA U^ 

+ 25m 2 L 2 ,U ( ^ - 25m 2 U (6E ^ - 2(5m 2 ) 2 L 2 *n (2 ^ + {5m 2 ) 2 Yi {Aa ^ + 2B 2 L 4b>1 U^ 
+ 2B 2 L 2 U^ - 2B 2 U^ - 2(5 2 ) 2 L 2 n^ + (E 2 ) 2 n( 4a ) - 2L 4bA U^ 
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n (p09) = n (p09) + 5m4aSm2ll (2**) + 5 m 4a B 2 U^ - 5m 4a U^ + 5m 2 B 4a U^ 

- 45m 2 B 2 L 2 U (2 * ) + 45m 2 L 2 n (4b ' 4 * ) - dm^ ^ - 2{5m 2 ) 2 L 2 U i2 ** ) + B 4a B 2 U (2) 

- B 4a U^ + AB 2 L 2 U^ - B 2 U^ - 2(B 2 ) 2 L 2 U^ - 2L 2 Tl^ 

n^ 0) = n (pl0) + 5m 4b 5m 2 U^ 2 ** ) + 5m 4b B 2 U^ - 5m 4b Tl^ b ' 4 * ) - 5m 2 Jm 2 B 2 U^ 

+ 5m 2 Jm 2 Ii i - AbM) - 5m 2 ^8m 2 ) 2 U^ + 5m 2 B 4b U (2 * ) - 5m 2 B 2 ,B 2 U i2) + 5m 2 B 2 ,U i4b) 
+ 25m 2 B 2 Tl^ + 5m 2 B 2 U^ - 25m 2 (B 2 ) 2 U^ - 5m 2 U^ c ^ - 5m 2 U^ B ^ 

- (5m 2 ) 2 B 2 *n^ - (5m 2 ) 2 B 2 U^ + (5m 2 ) 2 U (4b ' 2 *^ + B 4b B 2 U^ - B 4b U^ 

- B 2 U^ - fl 2 II< 6B > + 2(B 2 ) 2 U^ - (B 2 fnw 

n^ 1 ) = n( p11 ) - Sm 6F U^ + 25m 4a L 2 Il^ + 25m 2 L 4aA Il^ - 35m 2 (L 2 ) 2 U^ 

- B 6F TlW + 2B 4a L 2 Ii^ + 2B 2 L 4a ^ - W 2 (L 2 ) 2 Yi^ - 2L 4aA U.^ 

-2L 2 n( 6D ) + 3(L 2 ) 2 n( 4fe ) 

n^ 42 ) = n( p12 ) - 5m 6A Tl^ + 25m 4b(u) 5m 2 Tl^ + 25m 4b B 2 Yi (2 ^ - 5m 2 ^(5m 2 ) 2 Tl^ 

- 25m 2 Jm 2 B 2 U i2 * ) + 25m 2 B 4b{u) U {2) - 25m 2 B 2 ,B 2 U {2) + 25m 2 5 2 II (4b ' 2 * ) - 5m 2 (5 2 ) 2 n (2 * ) 

- 25m 2 U^ - (5m 2 ) 2 B 2 ^ + (5m 2 ) 2 U^ 2 * 2 ^ - B 6A U^ + 2B 4b B 2 U^ 

- 2B 2 U^ + (B 2 ) 2 U^ - (B 2 ) 3 U^ 

n (pi3) = n ( P i3) + 25m2LAa 2(u)Y[ m _ Sm 2 U^ + 2B 2 L 4a ^ - B 2 U^ - 2L 6DA U^ 

n (pi4) = n(p i4) + 25m 2 L 4bXu) TlM - 25m 2 L 2 *L 2 TlM + 25m 2 L 2 n( 4a ' 1 *) - 5m 2 (L 2 ) 2 n( 2 *) 

- 5m 2 n( 6G ' 2 *) + 2£? 2 L 4b , 2 n( 2 ) + 2fi 2 L 2 n( 4fl ) - 35 2 (L 2 ) 2 n( 2 ) - B 2 U^ 

- 2L eA2 U^ + 2L 4bA L 2 U^ - 2L 2 U^ + (L 2 ) 2 U^ 

n (pi5) = n ( P i5) _ 6mmY[ {2*) + 25m 2 L 4a , 2 U^ - B 6H ILW + 2B 2 L 4a ^ - 2L 4a>2 U^ 

ng6) = n(P i8) _ 2L 6HA nw + 2 L 4a , 2 L 2 n( 2 ) - L 4a , 2 n( 4a ) 
n (pi 7 ) = n(p i7) _ 2L 6 H, 2 n( 2 ) 
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n (pl8) = n ( P 18) _ 2L6G 2U (2) 



n (pi 9 ) = n(P i9) _ 2L 6Gi5 n( 2 ) + 2L 4b , 2 L 2 n( 2 ) - L 4M n^ + 2L 4(ljl L 2 n( 2 ) 

- L 2 n( 6F ) + (L 2 ) 2 n( 4fl ) - 2(L 2 ) 3 n( 2 ) 

n^ 2 °) = n( p2 °) - <Wn< 2 *) + 25m 4fe L 2 n( 2 *) + 5m 4a 5m 2 Jl^ + 8m 4a B 2 ,U^ 

- 5m 4 Jl^ 2 * ] - 25m 2 Jm 2 L 2 U {2 * ) + 5m 2 B 4a U {2 * ] - 25m 2 B 2 ,L 2 U i2) - 25m 2J B 2 L 2 II (2 * ) 
+ 25m 2 L 2 n( 4b ' 2 *) - B 6C U^ + 2B 4b L 2 TlW + B 4a B 2 Yi^ - B 4a U^ 

+ 2B 2 L 2 U^ - 2(B 2 ) 2 L 2 U^ - 2L 2 Tl^ 

= - 5m 6E U^ + 5m 4a{2if) 5m 2 U^ + 5m 4a B 2 U (2 ^ + 8m 2 B 4a(2 ^ 

- 25m 2 B 2 L 2 Jl {2) - 25m 2J B 2 L 2 n {2 * ) + 25m 2 L 4feil n {2 * ) + 25m 2 L 2 *II (46) - Sm^ 60 ' 3 ^ 

- 2(5m 2 )%*n( 2 *) - B 6E Il^ + B 4a B 2 Yl^ + 2B 2 L 4bA U^ + 2B 2 L 2 U^ 

- B 2 U^ - 2(B 2 ) 2 L 2 U^ - 2L 4btl U^ 

= U (p22 ^ - 5m 6B n {2 * } + 5m 4b{2if) 5m 2 U (2 ^ + 5m 4b 5m 2if U^ + 5m 4h B 2 Jl^ 

+ Sm^IlW - 5m 4b Tl^ b ' 2 * ) - 5m 2 Jm 2 B 2 *Tl {2) - 5m 2 Jm 2 B 2 U^ + 5m 2 Jm 2 Tl^ b ' 2 * ) 

- (8m 2 ,) 2 8m 2 U^ + 5m 2 B 4b{2)f) Yi {2) + 5m 2 B 4b U^ - 25m 2 B 2 ,B 2 U {2) + 5m 2 B 2 Ji {ib) 
+ 5m 2J B 2 n( 46 ' 2 *) - 5m 2 (i? 2 ) 2 n( 2 *) - 5m 2 U^ c ^ - (5m 2 ) 2 B 2 ^ - B 6B U^ 

+ 2B 4b B 2 U^ - B 4b U^ - B 2 ILW + (B 2 ) 2 U^ - (B 2 ) 3 U^ 

ngs) = n(P 23) _ 2L 6Gi3 n( 2 ) 

Iff > = n( p25 ) - 25m 2J B 2 L 2 ,n( 2 ) - 25m 2J B 2 L 2 n( 2 *) + 25m 2 B 2 Yi^ a ^ + 8m 2 L 4bm U^ 

+ 5m 2 L 4bM2 ^ + 25m 2 L 2 U^ - 5m 2 U^ E ^ - 5m 2 Yi^ E ^ - (5m 2 ) 2 L 2(UU) U^ 

- (5m 2 ) 2 L 2 n (2 ** ) + (5m 2 ) 2 n (4a ' 1 * 1 * ) + 25 2 L 4fejl n (2) + 25 2 L 2 n (4b) - 2B 2 Yi^ E) 

- 2{B 2 ) 2 L 2 ii^ + (5 2 ) 2 n( 4a ) - L 6A)1 n( 2 ) - L 2 U^ 
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n (p26) = n ( P 26) + Sm4aSm2U ^*) + Sm Aa B 2 U {2 ^ - 5m 4a U^ + 5m 2 B 4a U {2 ^ 

- 45m 2 B 2 L 2 U (2 * ) + 45m 2 L 2 n (4b ' 1 * ) - 5m 2 U (6D ' u) - 2{5m 2 ) 2 L 2 U i2 ** ) + B 4a B 2 U (2) 

- B 4a U^ + AB 2 L 2 U^ - B 2 Tl^ - 2(B 2 ) 2 L 2 Il^ - 2L 2 U^ 

n(ff ) = n( p2? ) + 5m 2 L 4 a,i(2*) n(2) + Sm 2 L 4ayl U^ - 5m 2 L 2 ,L 2 U^ + 5m 2 L 2 U^ a ^ 

- 5m 2 (L 2 ) 2 U^ - 8m 2 rt 6F ^ + 2B 2 L 4a ^ + B 2 L 2 U^ - 2B 2 (L 2 ) 2 U^ 

- B 2 n^ - l 6D ^ + L 4M L 2 n( 2 ) - L 4ail nw - l 2 u^ 
+ (L 2 ) 2 n^ 

n£f } = n (p28) + 5m 46 5m 2 n( 2 **) + 5m 4b B 2 U^ - 5m 4b Ii {ib ' u) - 5m 2 Jm 2 B 2 U^ 

+ 5m 2 Jm 2 rt 4b ^ - 5m 2 *(5m 2 ) 2 Tl^ + 5m 2 B 4b Tl^ - 5m 2 B 2 ,B 2 U^ + 5m 2 B 2 Jl^ 
+ 5m 2J B 2 n( 4fe ' 2 *) + 25m 2J B 2 n( 46 ' 1 *) - 25m 2 ( J B 2 ) 2 n( 2 *) - 5m 2 U^ c ^ - <5m 2 n( 6A ' 2 *) 

- (5m 2 ) 2 5 2 ,n( 2 *) - (5m 2 ) 2 5 2 n( 2 **) + (5m 2 ) 2 U^ b ' U2 ^ + 5 46 5 2 n( 2 ) - B 4b U^ 

- B 2 n^ - b 2 u^ + 2(B 2 ) 2 n(*» - (5 2 ) 3 n( 2 ) 

n (p29) = n ( P 2 9 ) + S m 2 L 4a , 2{2 ^ + 5m 2 L 4a , 2 U^ - 5m 2 U^ + 25 2 L 4a , 2 n( 2 ) 

- b 2 u^ - L 6E , 2 nW - L 4a , 2 nw 

n^°> = n( p3 °) + 5m 2 L 4bm U^ + 5m 2 L 45 , 2 n( 2 *) - 25m 2 L 2 ,L 2 n< 2 > + 5m 2 L 2 ,U^ 
+ 5m 2 L 2 rt 4a ^ - 5m 2 (L 2 ) 2 U^ - 5m 2 U {QG ^ + 2 J B 2 L 4b , 2 n( 2 ) + 2S 2 L 2 n^ 

- 35 2 (L 2 ) 2 n( 2 ) - b 2 u^ - L 6B)2 n( 2 ) - L 4b)2 U^ + 2L 4M L 2 n( 2 ) 

- L 4btl u^ - L 2 n< 6i5 > + (L 2 ) 2 nw 

n (p3i) = n(P 3i) + 5m 2 L 4a)1(1 ,)n( 2 ) + 5m 2 L 4ail n( 2 *) - 5m 2 L 2 *L 2 n( 2 ) + ^m^II^ 1 *) 

- 5m 2 (L 2 ) 2 n( 2 *) - 5m 2 n( 6F ' 5 *) + 2 J B 2 L 4a , 1 n( 2 ) + B 2 L 2 U^ - 25 2 (L 2 ) 2 n( 2 ) 

- b 2 u^ - L 6D , 3 n^ + L 4M L 2 n( 2 ) - L 4aA u^ - l 2 u^ 
+ (L 2 ) 2 iK 4fe ) 

n (p32) = n ( P 32) + 5mAaL2ifIl (2) + 6m4aL2ll (2*) _ ^n^' 1 *) - 25m 2 L 2 *L 2 n( 2 ) 

+ 25m 2 L 2 U^ a ^ - 25m 2 (L 2 ) 2 U^ + 2B 4a L 2 U^ - B 4a Il^ + 2B 2 L 2 U^ 

- AB 2 {L 2 ) 2 Yi^ - L 6C)1 n( 2 ) + 2L 4M L 2 n( 2 ) - 2L 2 U^ - L 2 U^ 

+ 2(L 2 ) 2 n^ 
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n (p33) = n ( P 33) _ L6G iU (2) _ ^^2) + AL4a lL2U m _ L4a lU (4a) 

- L 2 U^ + (L 2 ) 2 U^ - 2(L 2 ) 3 n( 2 ) 

n (p34) = n(p3 4) _ s m( , Dl i(2*) + 5 m46 L 2 n( 2 *) + 5m Mu) 5m 2 U (2 ^ + 5m 4aJ B 2 n( 2 *) 

- 5m 2 ^m 2 L 2 n( 2 *) + 5m 2 B Mu) U {2) - 5m 2 B 2 ,L 2 U^ - 5m 2 B 2 L 2 ,II (2) - 25m 2 fi 2 L 2 n (2 * ) 
+ 5m 2 L 4b ^ + 5m 2 L 2 Ji^ + <5m 2 L 2 n( 4fe ' 2 *) - 5m 2 U^ D ^ - (5m 2 ) 2 L 2 ^ 

- 5 6D n (2) + B 4b L 2 Yi {2) + 5 4a 5 2 n (2) + B 2 L 4bA U (2) + 25 2 L 2 n (4b) 

- B 2 n( 6D ) - 2(b 2 ) 2 l 2 tiw - L 4M n( 4fe ) - L 2 U^ 

n (p35) = n ( P 35) _ L6G 4U (2) _ ^^(2) + 2L4a>2L2 n(2) _ ^H) 

n (p36) = n ( P 36) _ L6Di2n ( 2 ) _ L6C2n (2) + 2L 4b;2 L 2 n( 2 ) + 3L 4ail L 2 n( 2 ) 

- L 4ajl n( 4a ) - L 2 n( 6G ) - L 2 n( 6F ) + 2(L 2 ) 2 n( 4a ) - 3(L 2 ) 3 n( 2 ) 

n r ^f) = + SmuLzJl® + 5m 4b L 2 U^ - 5m 4b U^ a ^ - 5m 2 Jm 2 L 2 ^ 

- 5m 2 Jm 2 L 2 U^ + 5m 2 Jm 2 rt 4a ^ - 25m 2 5 2 ,L 2 n( 2 ) + 5m 2 B 2 Ji^ - 5m 2 B 2 L 2 Ji^ 

- 5m 2 B 2 L 2 U^ + 5m 2 B 2 rt 4a ^ + 8m 2 L 4bm U^ + 5m 2 L 2 U {4b ^ - 5m 2 U^ E ^ 

+ 2B 4b L 2 n^ - B 4b n^ + B 2 L 4b ^ + 5 2 L 2 n( 4b ) - B 2 n^ 

- 2(B 2 ) 2 L 2 Il^ + (B 2 ) 2 U^ - L 6BA U^ - L 2 U^ 

nff) = II^ 38 ) + Sm 2 L 4aM4 ^ + <5m 2 L 4ail(3 *)n( 2 ) - 25m 2 L 2 *L 2 n( 2 ) + 5m 2 L 2 Ji^ 

- 5m 2 rt 6F ^ + 2B 2 L 4a ^ + B 2 L 2 U^ - 2B 2 (L 2 ) 2 U^ - B 2 U^ 

- L 6EA U^ - L eD , 5 n^ + 2L 4M L 2 n( 2 ) - L 4bA U^ 

n (p39) = n ( P 39) _ 5m6G n( 2 *) + 5m 4a L 2 U^ + 5m 2 L^ 2 U^ + 5m 2 L 4ajl n( 2 *) 

- 25m 2 (L 2 ) 2 n (2 * ) - £ 6G n (2) + 5 4a L 2 n (2) + 5 2 L 4M n (2) + 5 2 L 4ajl n (2) 

- 25 2 (L 2 ) 2 n( 2 ) - L 4b , 2 u^ - L 4aA nW _ i^um + 2 (L 2 ) 2 n( 4fe ) 
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Appendix C: Intermediate renormalization by if-operation 



This Appendix describes how UV-divergent subdiagrams obtained by applying the K- 
operation on the original unrenormalized vacuum-polarization functions are separated out. 

L^ v , -B^ v , etc. denote UV-divergent parts of renormalization constants L 2 , B 2 , etc. de- 
fined by the A-operation. Quantities with A attached in front are finite parts of the quan- 
tities. Note that A5m 2 = so that 5m 2 Y can be replaced by 5m 2 . Derivative amplitudes 



25] are denoted as L' 2 , B' 2 , etc. 



1. Fourth-order vacuum-polarization 



n (M =An (4a) +2L UV n (2) 

n (46) = An (4 fe ) + B UV n (2) + n (2*) 

2. Sixth-order vacuum-polarization 



n (6A) = An (6A) + 2 B VV n (46) _ {B UV )2 n (2) + 2 g^ n (46,l*) _ 2 g^ B UV n (2.) 

- (5m 2 ) 2 n( 2 **) 

n (6B) = An (6B) + 2 B VV n (46) _ (£? UV )2 n (2) + 2 ^ n (4 M *) _ 3 ^ 3UV n (2*) 

- (5m 2 ) 2 n( 2 **) 

n (6C) = An (6C) + B UV n (4 b ) _ S TJV ^/UV n (2) + 3UV n (2) _ ^/UV £ UV n (2*) 

+ 5m 2 n( 4b ' 2 *) - 5m 2 5m 2 u / n^ 2 *) + 5m u b v n( 2 *) 

U m = AU (6D) + B VV n (2) + 2 L UV n (46) _ 3 £ UV £ UV n (2) _ 3 ^ £ UV n (2*) 
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n (6£) = AU (6E) + B VW n (4a) _ ^UV £ UV n (2) + L UV n (46) _ ^UV £ UV n (2) 

+ l™l n (2) + 5m 2 n^' 1 *) - 5m 2 L 2 UV n( 2 *) 

n (6F) = An (6F) + L UV n (4a) _ 3 (L UV)2 n (2) + 3 ^UV n (2) 
n (6G) = An (6G) + 2 L UV n (4a) _ 3 (L UV )2 n (2) + 3 ^UV n (2) 

n(6^) = A n(6^) + 2Lu v 2 n( 2 ) 

3. Eighth-order vacuum-polarization 



n(p01) = An (p01) + 2 L UV n (6F) _ (L UV)2 n (4a) + 3 ( L UV)3 n (2) _ 4 ^UV L UV n (2) 
n (p02) =An ( P 02) +2jL UV 3n (2) 

n ( P 03) = An ( P 03) + 2 L UV n (6G) _ 3 (L UV )2 n (4a) + 4 (L UV )3 n (2) + 3 ^UV n (4a) 

-6LUV L uvn( 2 ) + 2^V 3 n( 2 ) 

n (p04) = An ( P 04) + 3 B VW n (6A) _ 3 (j g UV )2 n (46) + (j£ > UV )3 n (2) + ^ n (6A,3*) 

+ 2 5m 2 n^' 1 *) - 6 <5m 2 5 2 uv n^' 1 *) + 3 <5m 2 (5 2 uv ) 2 n< 2 *> - (5m 2 ) 2 n^ 1 * 3 *) 
- 2 {5m 2 f n^' 1 * 1 *) + 3 (<5m 2 ) 2 5 2 uv II^ 2 **) + {5m 2 f 

n ( P 05) = An ( P 05) + 2 S UV n (6B) + B VW n (6A) _ 3 (£? UV )2 n (46) + (j g> UV )3 n (2) 

+ 2 5m 2 n^' 1 *) + 5m 2 U^ 6A ^ - 4 5m 2 B 2 VV - 2 5m 2 B 2 VY 

+ 3 5m 2 (5 2 uv ) 2 n^ 2 *) - (5m 2 ) 2 n( 4M * 4 *) - 2 (<5m 2 ) 2 n^' 1 * 4 *) + 2 (<5m 2 ) 2 5 2 uv n^ 2 **) 

+ (5m 2 ) 2 B 2 VY n (2 ** } + (5m 2 f n (2 *** } 
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n ( P 06) = An ( P 06) + 2 £UV U (6E) _ (j g UV )2 n (4a) _ ^UV £ UV n (46) + ^UV (j g> UV )2 jj 

- V™ 5 2 uv + l'™ (5 2 uv ) 2 n< 2 ) + 2 l™ n( 4fe ) - 2 b™ n^ 2 ) 

+ 2 5m 2 II^' 1 *) - 2 5m 2 5 2 uv II^' 1 *) + 5m 2 L'™ E 2 UV 

+ 5m 2 L'™ 5 2 uv n< 2 *) - 2 5m 2 - (5m 2 ) 2 U^' 1 * 2 *) 

U (P07) = AU ( P 07) + 2 S UV U (6E) _ (j gUV )2 n (4a) + L »UV (£? UV )2 n (2) + ^UV n (6B) 

- 2 l 2 uv 5 2 uv n( 4b ) + l 2 uv (5 2 uv ) 2 n^ 2 ) - 2 Lyv (i/) b uv n(2) 

+ L™ 3 n^ 2 ) + 2 <5m 2 n( 6£ ' 5 *) - 2 5m 2 B™ - 2 5m 2 L™ n^ 4 *) 

+ 2 5m 2 L™ £ 2 UV n< 2 *> - (5m 2 ) 2 n^' 1 * 4 *) + (5m 2 ) 2 L 2 UV II^ 2 **) 

n ( P 08) = An ( P 08) + 2 S UV n (6B) _ (j gUV )2 n (4a) _ 3 ^/UV ^UV n (4*) + 3 ^UV ^UV) 

+ 2 n<*> - 2 5 2 UV n^ 2 ) + 2 5m 2 II^ 6 *) - 2 5m 2 B™ 
+ 2 5m 2 L'™ 5 2 uv n^ 2 *) - 2 5m 2 n( 2 *) - (5m 2 ) 2 n^' 1 * 3 *) 

n (p09) = An ( P 09) + B UV n (6D) + £UV n (46) _ ^UV £ UV n (2) + 3 £ UV n (6B) 

- 4 L 2 UV 5 2 uv Yi^ + 2 L 2 UV (i? 2 uv ) 2 11™ + 5m 2 Yi (%D ^ - 5m 2 B™ Ti^ 

- 4 5m 2 L™ n( 4M *) + 4 5m 2 L™ £ 2 UV + 2 (5m 2 ) 2 L 2 UV II^ 2 **) 
+ 5m™ - 5m™ B™ - 5m™ 5m 2 

n (plO) = An (plO) + B UV n (6C) + ^UV n (6B) _ B VV ^/UV n (4 b ) _ (£? UV )2 n (46) 

+ (5 2 uv ) 2 B'™ n^ 2 ) + B™ - B™ B™ U™ - 5m'™ B™ Il( 4fe > 4 *) 
+ 5m'™ (B™) 2 n^ 2 *) + 5m 2 U^ 6C ^ + 5m 2 U^ 6B ^ - 5m 2 B™ 

- 5m 2 B™ n^ 2 *) + 5m 2 B™ B'™ n< 2 *> - 5m 2 B™ Ii™ ] 

- 5m 2 5m™ II( 46 ' 4 *) + 5m 2 5m™ B™ + 5m 2 5m'™ B™ TL™^ 

- (5m 2 ) 2 Il( 4fe > 2 * 4 *) + (5m 2 ) 2 5m™ Tl^ + 5m™ - 5m™ B™ 
-5m™5m 2 U^ 

n (pll) = An (pH) + S UV n (2) + 2 L UV n (6D) _ 2 L UV S UV n (2) _ 3 (L UV)2 n (46) 

+ 3 (L 2 UV ) 2 B™ Yi^ + 2 L^i - 2 5 2 uv + 3 5m 2 {L™) 2 Yi^ 

- 2 5m 2 L™ A II( 2 *) - 2 5mY a v L 2 UV n^ 2 *) + 5m™ 
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n (pl2) = An (pl2) + 2 B UV n (6C) _ (j g UV )2 n (46) + (j gUV )2 ^UV n (2) _ fl W £ UV n (2) 

- s 4 u 6 y i;) 5 2 uv + 5 6 U A V n( 2 ) + W 2 /uv (£ 2 UV ) 2 n( 2 *) + 2 5m 2 n( 6C » 

- 2 5m 2 5 2 uv Il( 4fe » + 5m 2 5m™ B™ iK 2 *) + 5m 2 5m™ B™ Tl^ 

- {5m 2 f H^ 2 * 2 *) + (5m 2 ) 2 fcn™ n( 2 *) - 2 5m™ {u) 5m 2 - 5m™, } B ™ 
-5m™ v) B™Tl^+5m™Tl^ 

n (pl3) = An (pl3) + ^ UV n (6H) _ L UV (4#) s UV n (2) _ ^UV ^ £ UV n (2) 

+ 2LUV 4 n( 2 ) + ( 5m 2 n( 6H » 

n (pl4) = An (pl4) + B UV n (6G) + 2 L UV n (6E) _ 3 £ UV £ UV n (4a) + ^ UV jjXJW R UV n (2) 

+ l 2 uv l' 2 uv 5 2 uv n^ 2 ) - (l 2 uv ) 2 n( 4b ) + (l 2 uv ) 2 b™ 

rUV r UVtt(2) rUV R UV n (2) orUV rUVn(2) 

- L 46,2(4') B 2 11 - L 46,2(l') ^2 11 ~ 1 L 4b,l L 2 11 

+ 2 H ( 2 ) + 5m 2 n( 6G ' 2 *) - 2 5m 2 L™ n^' 1 *) + 5m 2 (L 2 UV ) 2 n< 2 *> 

n (pl5) = An (pl5) + B UV n (2) + 2 L UV n (4 6 ) _ 2 L UV 3 UV n (2) _ 3 ^ ^UV n (2.) 
n (pl6) = An (pl6) + L UV n (4a) _ 2 L UV ^ UV n (2) + 3 ^ n (2) 

n (pl7) = A n(p 17 )+2L™ 2n (2) 

n( P 18) =An ( P 18) +2L UV 2n (2) 

n (pl9) = An (pl9) + L UV n (6F) _ (L UV )2 n (4a) + 3 (L UV )3 n (2) + ^UV n (4a) 

- 2 L™ L™ n< 2 > - 2 L™ L™ n^ 2 ) + 2 L™ 5 n< 2 > 

n ( P 20) = An ( P 20) + S UV n (4 6 ) _ S UV ^/UV n (2) + 3UV n (2) + 3 L UV n (6C) 

- 2 L 2 UV 5 2 uv n(*> + 2 L 2 UV 5 2 uv 5 2 uv n^ 2 ) - 2 L 2 UV B™ 

- 5m™ B™ + 2 5m™ L™ B™ 11^ - 2 5m 2 L™ n^ 46 ' 2 *) 
+ 2 5m 2 5m™ L™ n^ 2 *) - 2 5m™ L™ Yi^ + 5m 4 u a v II^ 2 *) 

- 5m™ 5m™ Yi^ + 5m™ 
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n (p21) = An (p21) + £ UV n (6D) _ S UV #) s UV n (2) + 3UV n (2) _ ^UV £ UV n (46) 

+ l>™ (£ 2 UV ) 2 n( 2 ) - l'™ 5 2 uv + l' 2 uv (5 2 uv ) 2 n( 2 ) 

+ 2 - 2 S 2 uv + 5m 2 U^ D '^ + 5m 2 L'™ B™ Il( 2 *) 

+ 5m 2 L>™ B™ n^ 2 *) - 2 5m 2 L 4 U ^ n ( 2 *) - 5m™ 2 ^ 5m 2 Tl^ 
-5m™ 2l) B™U^+5m™U^ 

U (P22) = An ( P 22) + S UV n (6C) _ ^UV £/UV n (4 ) + ^ UV (j g,UV )2 n (2) 

- s 4 u b y 2 , } b™ n^ 2 ) + s 4 u 6 v - b™ b'™ iK 2 ) + b™ n ( 2 ) 

- <W 2 UV 5 2 uv n( 4b < 2 *) + 5m'™ B™ B'™ Il( 2 *) - 5m'™ B™ iK 2 *) 
+ 5m'™ 5m 2 u 7 5 2 uv Il( 2 *) + 5m 2 Tl^ - 5m 2 5m™ 

+ 5m 2 {5m™ f II( 2 *) - 5m™^ 5m 2 Tl^ - 5m™ 2l) B™ 
+ 5m™ n( 4fe » - 5m™ 5m™ + 5m™ 

n( P 23) =An ( P 23) +2L UV 3n (2) 

U (P24) = An ( P 24) + L UV n (6H) + ^UV n (4a) _ 3 ^UV ^ UV n (2) + ^ n (2) 

n (p25) = An ( P 25) + 2 S UV n (6B) _ (j gUV )2 n (4a) + ^//UV (£? UV )2 n (2) + ^UV n (6A) 

- 2 l 2 uv 5 2 uv n(*> + l™ (B™) 2 - l™ {4>) b™ n( 2 ) 

- ^1(2') B™ n( 2 ) + L™, n( 2 ) + 5m 2 U^ E ^ + 5m 2 U^ E ^ 

- 2 5m 2 £ 2 UV n^' 1 *) - 2 5m 2 L 2 UV II^ 1 *) + 2 5m 2 L 2 UV B™ 

- {5m 2 f n^' 1 * 1 *) + (5m 2 ) 2 L 2 uv Il( 2 **) 

n ( P 26) = An ( P 26) + B UV n (6D) + 3UV n (4 6 ) _ ^UV 3 UV n (2) + 3 L UV n (6A) 

- 4 L 2 UV 5 2 uv + 2 L 2 UV (B 2 UV ) 2 11™ + 5m 2 Yi {%D ^ - 5m 2 B™ Yi^ 

- 4 5m 2 L™ n^' 1 *) + 4 5m 2 L™ B™ li^ + 2 (5m 2 ) 2 L 2 UV n^ 2 ") 
+ 5m™ n^' 1 *) - 5m 4 u a v 5 2 uv n^ 2 *) - 5m™ 5m 2 
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n ( P 27) = An ( P 27) + B UV n (6F) + ^ UV n (6E) _ L UV £ UV n (4a) + ^ UV ^UV £ UV n (2) 

- (L 2 UV ) 2 n( 46 ) + (L 2 UV ) 2 b™ n( 2 ) - l™ n( 2 ) - l™^ 5 2 uv n( 2 ) 
+ l™ a n( 4b ) - l™, 5 2 uv n^ 2 ) + n^ 2 ) + 5m 2 n^ 1 *) 

- 5m 2 L™ n^' 1 *) + 5m 2 (L 2 UV ) 2 n^ 2 *) - 5m 2 L™, 11^ 

U (P28) = An ( P 28) + S UV n (6C) + ^UV n (6A) _ £ UV ^/UV n (46) _ (j gUV )2 n (4 6 ) 

+ (s 2 uv ) 2 fi 2 uv n< 2 ) + b™ _ B ™ b™ - 5m'™ b™ n^- 1 *) 

+ <W 2 UV (_B 2 UV ) 2 n^ 2 *) + 5m 2 Il^ c ^ + 5m 2 n( 6A » - 5m 2 B™ 

- 5m 2 B™ n^' 1 *) + 5m 2 B™ B 2 UV n( 2 *) - 5m 2 B™ 

- 5m 2 5m™ II^' 1 *) + 5m 2 5m™ B™ n^ 2 *) + 5m 2 5m'™ B™ 

- (5m 2 ) 2 nW* 2 *) + (5m 2 ) 2 5m™ + 5m™ n^ 1 *) - 5m™ B™ 
-5m™5m 2 U^ 

n ( P 29) = An ( P 29) + B UV n (6if) _ L UV ^ £ UV n (2) + ^UV n (46) _ ^UV £ UV n (2) 

+ L™ 2 n( 2 ) + 5m 2 n^ 1 *) - 5m 2 L 4 u v n ( 2 *) 

n ( P 30) = An ( P 30) + B UV n (6G) _ L ,UV £ UV n (4a) + ^ UV n (6J5) _ £ UV £ UV n (4a) 

+ 2 l 2 uv l 2 uv e 2 uv n( 2 ) - (L 2 UV ) 2 n( 4b ) + (L 2 UV ) 2 b™ 

- lZ (30 5 2 uv n< 2 > + n<"> - 5 2 uv n( 2 ) + l™ 

- 2 L™ L™ n( 2 ) + Lyv 2 n^ 2 ) + 5m 2 n^' 1 *) - 5m 2 L™ n^' 1 *) 
+ 5m 2 (L™) 2 n^ 2 *) - 5m 2 L™ 

n (p31) = An (p31) + £ UV n (6F) + L UV U (6E) _ ^ UV £ UV n (4a) + L UV ^/UV £ UV n (2) 

- (L 2 UV ) 2 n(*> + (l 2 uv ) 2 b™ - l™ n^ 2 ) - L u a v 1(10 5 2 uv n^ 2 ) 

+ L™ n«*> - L™ B™ n< 2 > + n^ 2 ) + 5m 2 

- 5m 2 L 2 UV n^' 1 *) + 5m 2 (L 2 UV ) 2 n^ 2 *) - 5m 2 L™, 

n (p32) = An ( P 32) + B TJV n (4a) _ L /UV ^UV n (2) + 3 L UV U (6E) + ^ UV n (6D) 

- 2 l 2 uv 5 2 uv n( 4fl ) - l 2 uv 5 4 u a v n^ 2 ) + 2 l 2 uv l' 2 uv 5 2 uv n^ 2 ) 

- 2 (L 2 UV ) 2 + 2 (L 2 UV ) 2 B™ - 2 L™ L 2 UV n^ 2 ) + II^ 

- 2 5m 2 L 2 UV II^' 1 *) + 2 5m 2 (L 2 UV ) 2 Il( 2 *) + 5m™ Ii {Aa ^ - 5m™ L™ Yi^ 
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n ( P 33) = An ( P 33) + L UV n (6F) _ (L UV)2 n (4a) + 3 (L UV)3 n (2) + ^UV n (4a) 

_ 4 rUV r UV tt(2) , rUV tt(2) , rUV tt(2) 
4 ^4o,l ^2 11 + ^6G,1 11 + ^6F,1 11 

n (p34) = An ( P 34) + B UV n (6D) _ ^ £ UV n (2) + ^UV n (2) _ ^/UV £ UV n (4ft) 

+ 4 UV (£? 2 UV ) 2 n( 2 ) + l™ n( 6C ) - l 2 uv # 2 UV + l™ 5 2 uv 5 2 uv n( 2 ) 

- l™ s 4 u & v n( 2 ) + l™ n< 46 > - s 2 uv n^ 2 ) + 5m 2 uv l 2 uv i? 2 uv n^ 2 *) 

+ 5m 2 n( 6D ' 2 *) + 5m 2 L'™ 5 2 uv n^ 2 *) - 5m 2 L™ II^' 2 *) 

- 5m 2 L™ n^ 2 *) + 5m 2 5m 2 u / L 2 UV n^ 2 *) - 6m™ L™ n^ 2 *) 

- 5m 2 n( 2 *) - 5m™ m B ™ n( 2 *) + 5m™ 

n ( P 35) = An ( P 35) + L UV U (6H) _ 3 ^UV £ UV n (2) + ^ n (2) + ^UV n (2) 
n (p36) = An ( P 36) + L UV n (6G) + ^ UV n (6F) _ 3 (L UV )2 n (4a) + 3 ^ UV )3 n (2) 

- 2 l 2 uv n^ 2 ) + l™ n<*o - 3 l™ n( 2 ) + n ( 2 ) 
+ l™ 2 u^ 

n (p37) = An ( P 37) + B UV n (6£0 _ £ UV £/UV n (4a) + ^UV n (4a) + ^/UV £ UV £/UV n (2) 

- U™ s 4 u 6 v n^ 2 ) + l™ n( 6C ) - l 2 uv 5 2 uv + l 2 uv s 2 uv 5 2 uv n^ 2 ) 

- L™ 5,7 n< 2 > - L™ m B™ n< 2 > + n ( 2 > - 5m 2 uv 5 2 uv n^ 1 *) 
+ 5m'™ L™ £ 2 UV II( 2 *) + 5m 2 U^ E ^ - 5m 2 L™ II^ 2 *) 

- 5m 2 5m™ Tl^ + 5m 2 5m™ L™ Tl^ + 5m™ 
-5m™L™U^ 

U (P38) = An ( P 38) + B UV n (6F) _ ^UV 3 UV n (4a) + 3 ^ UV jjVW R UV n (2) 

+ n(-) - 2 L 2 UV n( 2 ) - L 4 u V (4/) 5 uv n (2) _ L w (y) s uv n (2) 
+ L™,U^ + L™ 5 Il^ + 5m 2 U^ 

n ( P 39) = An ( P 39) + B VV n (2) + L UV n (6D) _ ^ UV ^UV n (2) _ 3 (L UV )2 n (4 6 ) 

+ 2 (L 2 UV ) 2 5 2 uv n^ 2 ) + l™ _ L yv 5 2 uv n^ 2 ) + n( 4fe ) 

- L™, B™ li^ + 2 <5m 2 (L 2 UV ) 2 n^ 2 *) - 5m 2 L™ 

- 5m 2 L™ A - 5m™ L™ + 5m™ Tl^ 
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Appendix D: Divergence structure of renormalization constants of sixth and lower 
orders 

Throughout this Appendix L n , B n , and 8rh n denote quantities obtained by removing the 
overall UV divergences of L n , B n , and 5m n by the A"-operation. They may still have sub- 
diagram UV divergences which are subtracted by subdiagram if-operations. The resulting 
UV-finite quantities are denoted as L^, and 5m^. These quantities may have IR diver- 
gences, which are subtracted by .^-subtraction and /-subtraction. These operations create 
UV- and IR- finite quantities which are denoted as AL n , AB n , and ASm n . 

1. Second-order renormalization constants 



L 2 — L 2 + L 2 , L 2 — L 2 — J 2 , 

B 2 = B 2 V + B 2 , B^ — B 2 — —I 2 + AB 2 , 

ALB 2 = Lf + Bf = AB 2 , 

B 2 * — — 27^2*5 L 2 * — I 2 * -\- AL 2 * , 

B 2 ** = — 2(21/2(1*1*) +-^2(1*2*))) 

5m 2 * = 5m 2 y + I 2 + ASm 2 * . 
2. Fourth-order renormalization constants 



Sm Aa = 5m 4a 

5m^ b = 5iri4b — Sfh 2 , B 2 Y — 5m 2 * 5m 2 
A5m 4a = 5mf a 
A5m 4b = Sm^ b 
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A<5m 4 = A5m 4a + A5m 4b 



75R r> or UV r> 

-D 4a — &4a — * L 2 £>2 

B\\ = B Ab — B 2 , B 2 Y — 5m 2 B 2 * 

rR _ 7 7 rUV 

^4a,l _ ^4a,l _ ^2 ^2 

rR _ 7 

-^46,1 = ^46,1 — L 2 i B 2 Y — 5m 2 L 2 * 



rR _ r r rUV 

-^46,2 _ -^46,2 — L>2 i-> 2 



ALB 4a = Bf a + 2 Lf a l + L 4a2 
ALB 4b = Bf b + 2 L 4b;1 + L% 2 - L 2 Bf - (L 2 ) 2 
ALB 4 = ALB 4a + ALB 4fe 
3. Sixth-order renormalization constants 



^ m aA = $ m 6A - -B 2 v £ra 46 (i/) - -83 v 5m 4fe ( 3 /) - 2 <5m 2 5m 4b (i*) + Sm 2 „ (-B 2 v ) 2 + 8m 2if , 5m 2 B. 



?uv . 



?uv 



UVn2 



? uv 



+ 8m 2if , 5m 2 B 2 Y + 5m 2 ** (8m 2 



^ m 6B = 5m &B - B 2 V 5m 4b ( 2 ') - 5m 2 5m 4b(2 *) - 5m 2/ B^ b + 5m 2 , B^ V B^ V 
— 5m 2 * 5m 4b Y + 5m 2lf 5m 2 Y B 2 V + 5fh 2lf 5m 2 ? Sm 2 



5m^ c = 5m 6C - 2 L 2 V 5m 4b - 5m 2 , B^ + 2 5m 2 , L 2 y 5 2 uv - 5m 2if 5m v J 



+ 2 5m,2* Sm 2 L 2 Y 
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Smf D = SrriQ D - B^ v 5m Aa{v) 
+ 5m 2 * Sm 2 



L^ v 8m Ab - 5m 2 5m Aa ( U ) + Sm 2 , L^ v 



5mf E = 5m 6E - B^ v 5m Aa{2l) - Sm 2 8m Aa{2 * 



5mf F = 5rriQF — 2 L^ v 5m Aa 



5m% = 5m 6G - L 2 Y 5m Aa 



5m% H = 5m m 



A5m 6A = 5m^ A 



A5m 6B = 5mf B - L 2 5mf b 



A5m 6C = Smf c - L 2 5m Aa 



A5m 6D = 5m^ D 



A5m 6E = Sm^ E 



A5m 6F = Sm^ F 



A5m 6G = 5mf G 



A5m m = 5mf H 



ASm e = Sm^ A + 5m^ B + 5m^ c + 2 5m^ D + 8m% E + 5m^ F + 2 5m^ G + 5m^ H - L 2 8m^ a 
- L 2 5m^ b 



44 



Bqa — Bqa — B 2 y -046(1') — B 2 y -846(3') + B 2// (B 2 V ) 2 — 2 5m 2 £46(1*) + B^ v -B 2 *, 
+ 5m 2 5 2 uv 5 2 *, + (5m 2 ) 2 B 

2** 

= Bqb — -B 2 * Sm 4b v — B>2 Y B ib (2') — B 2 , B^ b + B 2 , B^ v B^ v 

- 5m 2 B Ab (2*) + Sm 2 ^ B 2 V B 2if + 5m 2 ^ 5m 2 B 2 * 

B^c = B%c — -B2* 5m u a V — B 2 , B^ — 2 L 2 V B Ab + 2 L 2 V B 2l B 2 y 
+ 2 5m 2 L^v £ 2 , 

R R _ 5 R UV 5 d rUV rUV 5 , rUV n R UV 

-°6D — — -D 2 -t>4a(l') — -°2 -^4b,l ~~ L 2 -°46 + J-"2 &2i -t> 2 

+ £ 2 5 2 uv - 5m 2 5 4a(H + 5m 2 Lf 1 B 2 * 

R R _ d R UV r 9 R 7~ uv i 7~ uv R R uv i r uv r r uv 

-°6E — -°6S — -D 2 -D4a(2') ~ ^ H 2 L 4bl + L 2l JJ 2 £> 2 + L 2 , t5 2 B 2 

— 5m 2 B 4a ^) 

Bqf = Bqf — 2 B 2 L\^ x — 2 L 2 V B Aa + 3 (£ 2 V ) 2 B 2 

RR _ 5 R rUV n rUV rUV n , 9 /rUV\2n 

-°6^ - -D6tf - ^ ^2 -^4a,2 

-^^4,1 = ^6A,1 — -B^ -^46,1((1')') — ^46,1(3') + L 2 „ (B 2 Y ) 2 — 5m 2 Z/46,l((l')*) — 5m 2 £46,1(3*) 

+ <5m 2 L 2 */ + 5m 2 L 2 */ + (Sm 2 ) 2 L 2 (i*i*) 

-^£L4,2 = -^6A,2 — -B^ ^46,2(3') — ^2^ ^46,1 + -^2/ -^ V -B^ — <^ m 2 ^46,2(1*) 

+ 5m 2 L 2 * 

-^6A,3 = L 6 a,s — 2 -B^ £46,1(1') + L 2 „ (B 2 Y ) 2 — 2 5m 2 £46,1(1*) + 2 ^2 £2*/ -B^ 
+ (5m 2 ) 2 L 2 (i*(i/ );(i ) 

45 



7-R _ t rUV t t Am UV j R UV , r R UV R UV 

-^65,1 ~~ -^65,1 — E>2 -^46,1(2') ~ -^2* 0m 4b — L21 n Ab + L21 D 2l &2 

- 5m 2 -^46,1(2*) + ^ m 2^ L 2* B 2 V + 5m 2 Y 5m 2 L 2 * 
^gb,2 = Lqb,2 — B 2 Y i-4&,2((2')') — L 2 + L 2 L 2 ^ B 2 Y — 5m 2 £46,2(2* 

^65,3 — ^65,3 - ^2 L *b,2 ~ ^2 -^ 4 fe,2 + L 2 {^2 ) 



tR t t A™ uv OfUVr t R UV , 9 r rUV nl 

-^6(7,1 — ^6C,1 ~ L 2* <jm Aa — I L 2 _L 4fe)1 — L 2l -D 4a + Z _L 2/ -L 2 £> 2 

UV r 
v 2 -^2* 



+ 2 <5m 2 L 2 



^6C,2 — Lqc,2 - L 2 Y L Ab>2 - L 2 L A ^ + L 2 (L 2 V ) 2 

T R _ 7 7 rUV 

-^60,3 - L GC,3 - ^2 -^4a,2 

^6D,1 = Lqd,1 ~ B 2 Y I/4o,l((l')') ~~ ^2^ ^46,1 + ^2/ -^ 2 V ~~ ^4a,l((l')* 

+ 5m 2 L 2 V L 2 * 

rR _7 rUV7 rUV7 ,7 /tUVn2 

-^6.0,2 — L §D,2 — L* 2 L A a : i — L 2 L Abj2 + Li 2 \L 2 ) 

^6D,3 = Lqd,3 — B 2 Y Li4a,\{V) ~ L 2 ^ Lib,l + L 2/ L 2 W B 2 V — 5m 2 -L4a,l(l*) 

+ 5m 2 L 2 V L 2 * 

^6D,4 = Lqd,4 — B 2 Y Lia,2{V) ~ $ m 2 ^4o,2(l*) 

L(5D,5 = -^6D,5 — ^4a,l(3') — ^2 ^4^1 + ^2 -^2^ — ^ m 2 ^4a,l(3*) 

-^6£,i = -^6E,i — B 2 V L Aa: i(2') — L 2 + L 2 L 2 Y B 2 Y — 5m 2 £40,1(2*) 

^6E,2 = L§E,2 — B 2 V L 40j2 ((2')') — ^2 -^4a,2(2*) 
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j-ti _ t r uv r 

-^6£,3 — ^6E,3 — L>2 ^4a,2 

rR _7 rUV? 7 rUV , 7 / r UVn2 

^6F,1 — ^6F,1 — ^2 ^4a,l _ -^2 -^4a,l + -^2 1-^2 J 

r R _ 7 r UV 7 

-^6F,2 — ^6F,2 — 1^2 -^4a,2 

rR _7 o rUV7 ,7 frUVN2 

-^6G,1 = ^6G,1 - ^4a,l - ^2 £4^1 + L 2 {L^Y 

^60,2 — L§G,2 

L 6G,3 = L 6G,3 

^6G,4 = Lqg,A — Lia,2 

rR _7 ruv7 7 rUV , 7 frUV\2 

-^60,5 — -^60,5 — ^2 ^4a,l _ ^2 -^46,2 + L 2 \^2 ) 

rR _ 7 7 rUV 

^611,2 — LqH,2 
L m,3 = L GH,3 

ALB 6A = B* A + 2 Lq AA + 2 L^ A 2 + Lq A 3 - 2ALB 2 L^ b l 

ALBqb = Bf B + 2 L^ B l + 2 Lg B2 + L^ B 3 - ALB Ab - ALB 2 L% 2 

AL£> 6 c = B& c + 2 + 2 Lg C 2 + L^ C 3 - ALB ia Lf 

ALBq D = Bf D + L^ D l + L^ D 2 + L^ D 3 + L^ DA + L 6D5 — ALB 2 Lf^ 
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ALBqe — Bf E + 2 Lq E>1 + 2 L^ E2 + L^ E3 - ALB 2 L^ a 2 



ALBqf — B^ F + 2 L^ F 1 + 2 2 + 



6F,3 



ALBq G — Bf G + Lg G1 + Lg G 2 + L 6G)3 + -^^ Gj4 + -^6G,£ 



ALBq H — Bf H + 2 + 2 L^ H2 + Lg H 3 



/3=A 

where A^ = \b = A G = A^ = \p = A# = 1, and A^ = A G = 2. 
AL 6 and AB 6 defined in Ref. [3| are related to ALBq through 



(Dl) 



ALB 6 = AL 6 + AB 6 + AL 4 AB 2 + A5m 4 B 2 * [/]. 



(D2) 
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